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ABSTRACT 



Cosmology now provides unambiguous, quantitative evidence for new 
particle physics. I discuss the implications of cosmology for supersym- 
metry and vice versa. Topics include: motivations for supersymmetry; su- 
persymmetry breaking; dark energy; freeze out and WIMPs; neutralino 
dark matter; cosmologically preferred regions of minimal supergravity; 
direct and indirect detection of neutralino s; the DAM A and HEAT sig- 
nals; inflation and reheating; gravitino dark matter; Big Bang nucleosyn- 
thesis; and the cosmic microwave background. I conclude with specula- 
tions about the prospects for a microscopic description of the dark universe, 
stressing the necessity of diverse experiments on both sides of the particle 
physics/cosmology interface. 
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1 Introduction 

Not long ago, particle physicists could often be heard bemoaning the lack of unam- 
biguous, quantitative evidence for physics beyond their standard model. Those days 
are gone. Although the standard model of particle physics remains one of the great 
triumphs of modern science, it now appears that it fails at even the most basic level — 
providing a reasonably complete catalog of the building blocks of our universe. 

Recent cosmological measurements have pinned down the amount of baryon, mat- 
ter, and dark energy in the universe.^ In units of the critical density, these energy 
densities are 

tt B = 0.044 ±0.004 (1) 

Onatter = 0.27 ± 0.04 (2) 

tt A = 0.73 ±0.04, (3) 
implying a non-baryonic dark matter component with 

0.094 < tt DM h 2 < 0.129 (95% CL) , (4) 

where h ~ 0.71 is the normalized Hubble expansion rate. Both the central values and 
uncertainties were nearly unthinkable even just a few years ago. These measurements 
are clear and surprisingly precise evidence that the known particles make up only a 
small fraction of the total energy density of the universe. Cosmology now provides 
overwhelming evidence for new particle physics. 
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At the same time, the microscopic properties of dark matter and dark energy are 
remarkably unconstrained by cosmological and astrophysical observations. Theoretical 
insights from particle physics are therefore required, both to suggest candidates for dark 
matter and dark energy and to identify experiments and observations that may confirm 
or exclude these speculations. 

Weak-scale supersymmetry is at present the most well-motivated framework for 
new particle physics. Its particle physics motivations are numerous and are reviewed in 
Sec.|2l More than that, it naturally provides dark matter candidates with approximately 
the right relic density. This fact provides a strong, fundamental, and completely inde- 
pendent motivation for supersymmetric theories. For these reasons, the implications of 
supersymmetry for cosmology, and vice versa, merit serious consideration. 

Many topics lie at the interface of particle physics and cosmology, and supersym- 
metry has something to say about nearly every one of them. Regrettably, spacetime 
constraints preclude detailed discussion of many of these topics. Although the discus- 
sion below will touch on a variety of subjects, it will focus on dark matter, where the 
connections between supersymmetry and cosmology are concrete and rich, the above- 
mentioned quantitative evidence is especially tantalizing, and the role of experiments 
is clear and promising. 

Weak-scale supersymmetry is briefly reviewed in Sec. |2] with a focus on aspects 
most relevant to astrophysics and cosmology. In Sees. |3] and |U the possible roles of 
neutralinos and gravitinos in the early universe are described. As will be seen, their 
cosmological and astrophysical implications are very different; together they illustrate 
the wealth of possibilities in supersymmetric cosmology. I conclude in Sec. |5] with 
speculations about the future prospects for a microscopic understanding of the dark 
universe. 

2 Supersymmetry Essentials 
2.1 A New Spacetime Symmetry 

Supersymmetry is an extension of the known spacetime symmetries. 3 The spacetime 
symmetries of rotations, boosts, and translations are generated by angular momentum 
operators L { , boost operators K { , and momentum operators P^, respectively. The L and 
K generators form Lorentz symmetry, and all 10 generators together form Poincare 
symmetry. Supersymmetry is the symmetry that results when these 10 generators are 
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further supplemented by fermionic operators Q a . It emerges naturally in string theory 
and, in a sense that may be made precise, 4 is the maximal possible extension of Poincare 
symmetry. 

If a symmetry exists in nature, acting on a physical state with any generator of 
the symmetry gives another physical state. For example, acting on an electron with a 
momentum operator produces another physical state, namely, an electron translated in 
space or time. Spacetime symmetries leave the quantum numbers of the state invariant 
— in this example, the initial and final states have the same mass, electric charge, etc. 

In an exactly supersymmetric world, then, acting on any physical state with the 
supersymmetry generator Q a produces another physical state. As with the other space- 
time generators, Q a does not change the mass, electric charge, and other quantum 
numbers of the physical state. In contrast to the Poincare generators, however, a su- 
persymmetric transformation changes bosons to fermions and vice versa. The basic 
prediction of supersymmetry is, then, that for every known particle there is another 
particle, its superpartner, with spin differing by |. 

One may show that no particle of the standard model is the superpartner of an- 
other. Supersymmetry therefore predicts a plethora of superpartners, none of which 
has been discovered. Mass degenerate superpartners cannot exist — they would have 
been discovered long ago — and so supersymmetry cannot be an exact symmetry. The 
only viable supersymmetric theories are therefore those with non-degenerate superpart- 
ners. This may be achieved by introducing supersymmetry-breaking contributions to 
superpartner masses to lift them beyond current search limits. At first sight, this would 
appear to be a drastic step that considerably detracts from the appeal of supersymmetry. 
It turns out, however, that the main virtues of supersymmetry are preserved even if such 
mass terms are introduced. In addition, the possibility of supersymmetric dark matter 
emerges naturally and beautifully in theories with broken supersymmetry. 

2.2 Supersymmetry and the Weak Scale 

Once supersymmetry is broken, the mass scale for superpartners is unconstrained. 
There is, however, a strong motivation for this scale to be the weak scale: the gauge 
hierarchy problem. In the standard model of particle physics, the classical mass of the 
Higgs boson (m\) receives quantum corrections. (See Fig. [TJ) Including quantum 
corrections from standard model fermions f L and f R , one finds that the physical Higgs 
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Fig. 1 . Contributions to the Higgs boson mass in the standard model and in supersym- 
metry. 



boson mass is 



1 



< = K)o - j^ XA + ■ ■ ■ . (5) 

where the last term is the leading quantum correction, with A the Higgs-fermion cou- 
pling. A is the ultraviolet cutoff of the loop integral, presumably some high scale well 
above the weak scale. If A is of the order of the Planck scale ~ 10 19 GeV, the classical 
Higgs mass and its quantum correction must cancel to an unbelievable 1 part in 10 34 to 
produce the required weak-scale m^. This unnatural fine-tuning is the gauge hierarchy 
problem. 

In the supersymmetric standard model, however, for every quantum correction with 
standard model fermions f L and f R in the loop, there are corresponding quantum cor- 
rections with superpartners fi and The physical Higgs mass then becomes 

ml = K)o-^A^ +I ^A 2 A 2 + ... 

~ {m 2 h )o + Y^( m / ~ m /) ln ( A /™/) > ( 6 ) 

where the terms quadratic in A cancel, leaving a term logarithmic in A as the leading 
contribution. In this case, the quantum corrections are reasonable even for very large 
A, and no fine-tuning is required. 

In the case of exact supersymmetry, where mj = mf, even the logarithmically di- 
vergent term vanishes. In fact, quantum corrections to masses vanish to all orders in 
perturbation theory, an example of powerful non-renormalization theorems in super- 
symmetry. From Eq. ©, however, we see that exact mass degeneracy is not required 
to solve the gauge hierarchy problem. What is required is that the dimensionless cou- 
plings A of standard model particles and their superpartners are identical, and that the 
superpartner masses be not too far above the weak scale (or else even the logarithmi- 
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Fig. 2. Neutral particles in the supersymmetric spectrum. Mi, M 2 , /i, and m 3 / 2 
are unknown weak-scale mass parameters. The Bino, Wino, and down- and up-type 
Higgsinos mix to form neutralinos. 

cally divergent term would be large compared to the weak scale, requiring another fine- 
tuned cancellation). This can be achieved simply by adding supersymmetry-breaking 
weak-scale masses for superpartners. In fact, other terms, such as some cubic scalar 
couplings, may also be added without re-introducing the fine-tuning. All such terms 
are called "soft," and the theory with weak-scale soft supersymmetry-breaking terms is 
"weak-scale supersymmetry." 

2.3 The Neutral Supersymmetric Spectrum 

Supersymmetric particles that are electrically neutral, and so promising dark matter 
candidates, are shown with their standard model partners in Fig. El In supersymmetric 
models, two Higgs doublets are required to give mass to all fermions. The two neutral 
Higgs bosons are H d and H u , which give mass to the down-type and up-type fermions, 
respectively, and each of these has a superpartner. Aside from this subtlety, the super- 
partner spectrum is exactly as one would expect. It consists of spin sneutrinos, one for 
each neutrino, the spin | gravitino, and the spin | Bino, neutral Wino, and down- and 
up-type Higgsinos. These states have masses determined (in part) by the corresponding 
mass parameters listed in the top row of Fig. El These parameters are unknown, but are 
presumably of the order of the weak scale, given the motivations described above. 
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Fig. 3. Proton decay mediated by squark. 



(7) 



The gravitino is a mass eigenstate with mass m 3 / 2 . The sneutrinos are also mass 
eigenstates, assuming flavor and i?-parity conservation. (See Sec. 12.41 ) The spin | 
states are differentiated only by their electroweak quantum numbers. After electroweak 
symmetry breaking, these gauge eigenstates therefore mix to form mass eigenstates. In 
the basis (—IB, — iW 3 , H d , H u ) the mixing matrix is 

/ Mi -M z cosps w M z sin[3s w \ 

M 2 M z cosf3c w -M z sm(3c w 

—Mz cos (3 sw Mz cos (3 cw —\i 

i M z sin j3 sw —M z sin (3 cw — yU J 

where cw = cos^vi/, Svf = sin^, and (3 is another unknown parameter defined 
by tan/5 = (H u )/(H d ), the ratio of the up-type to down-type Higgs scalar vacuum 
expectation values (vevs). The mass eigenstates are called neutralinos and denoted 
{x = Xii X2, X3; Xa}, in order of increasing mass. If Mi M 2 , |/i|, the lightest 
neutralino x nas a mass of approximately Mi and is nearly a pure Bino. However, for 
Mi ~ M 2 ~ x is a mixture with significant components of each gauge eigenstate. 

Finally, note that neutralinos are Majorana fermions; they are their own anti- 
particles. This fact has important consequences for neutralino dark matter, as will be 
discussed below. 



2.4 i?-Parity 

Weak-scale superpartners solve the gauge hierarchy problem through their virtual ef- 
fects. However, without additional structure, they also mediate baryon and lepton num- 
ber violation at unacceptable levels. For example, proton decay p — > n°e + may be 
mediated by a squark as shown in Fig. [3j 

An elegant way to forbid this decay is to impose the conservation of i?-parity 
R p = (— i) 3 (-B-L)+25^ w ] iere anc i g are baryon number, lepton number, and 
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spin, respectively. All standard model particles have R p = 1, and all superpartners 
have R p = — 1. i?-parity conservation implies TIR P = 1 at each vertex, and so both 
vertices in Fig.|3]are forbidden. Proton decay may be eliminated without i?-parity con- 
servation, for example, by forbidding B or L violation, but not both. However, in these 
cases, the non-vanishing i?-parity violating couplings are typically subject to stringent 
constraints from other processes, requiring some alternative explanation. 

An immediate consequence of i?-parity conservation is that the lightest supersym- 
metric particle (LSP) cannot decay to standard model particles and is therefore stable. 
Particle physics constraints therefore naturally suggest a symmetry that provides a new 
stable particle that may contribute significantly to the present energy density of the 
universe. 

2.5 Supersymmetry Breaking and Dark Energy 

Given imparity conservation, the identity of the LSP has great cosmological impor- 
tance. The gauge hierarchy problem is no help in identifying the LSP, as it may be 
solved with any superpartner masses, provided they are all of the order of the weak 
scale. What is required is an understanding of supersymmetry breaking, which governs 
the soft supersymmetry-breaking terms and the superpartner spectrum. 

The topic of supersymmetry breaking is technical and large. However, the most 
popular models have "hidden sector" supersymmetry breaking, and their essential fea- 
tures may be understood by analogy to electroweak symmetry breaking in the standard 
model. 

The interactions of the standard model may be divided into three sectors. (See 
Fig. |U) The electroweak symmetry breaking (EWSB) sector contains interactions in- 
volving only the Higgs boson (the Higgs potential); the observable sector contains in- 
teractions involving only what we might call the "observable fields," such as quarks q 
and leptons /; and the mediation sector contains all remaining interactions, which cou- 
ple the Higgs and observable fields (the Yukawa interactions). Electroweak symmetry 
is broken in the EWSB sector when the Higgs boson obtains a non-zero vev: h — > v. 
This is transmitted to the observable sector by the mediating interactions. The EWSB 
sector determines the overall scale of EWSB, but the interactions of the mediating sec- 
tor determine the detailed spectrum of the observed particles, as well as much of their 
phenomenology. 

Models with hidden sector supersymmetry breaking have a similar structure. They 
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Fig. 4. Sectors of interactions for electroweak symmetry breaking in the standard model 
and supersymmetry breaking in hidden sector supersymmetry breaking models. 



have a supersymmetry breaking sector, which contains interactions involving only 
fields Z that are not part of the standard model; an observable sector, which contains all 
interactions involving only standard model fields and their superpartners; and a media- 
tion sector, which contains all remaining interactions coupling fields Z to the standard 
model. Supersymmetry is broken in the supersymmetry breaking sector when one or 
more of the Z fields obtains a non-zero vev: Z — > F. This is then transmitted to the 
observable fields through the mediating interactions. In contrast to the case of EWSB, 
the supersymmetry-breaking vev F has mass dimension 2. (It is the vev of the auxiliary 
field of the Z supermultiplet.) 

In simple cases where only one non-zero F vev develops, the gravitino mass is 

where M, = (8wG N )' ,/2 ~ 2.4 x 10 ls GeV is the reduced Planck mass. The standard 
model superpartner masses are determined through the mediating interactions by terms 
such as ^ 

c ijTirfifj and c «F Afl ^' (9) 

where q,- and c a axe constants, f { and A a are superpartners of standard model fermions 
and gauge bosons, respectively, and M m is the mass scale of the mediating interactions. 
When Z — > F, these terms become mass terms for sfermions and gauginos. Assuming 
order one constants, 

F 

m 7Tt A — — — . (10) 
In supergravity models, the mediating interactions are gravitational, and so M m ~ 
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M*. We then have 

F 

m 3/2 ,m f ~,m x ~ — , (11) 

and \[F ~ x/Mweak-M* ~ 10 10 GeV. In such models with "high-scale" supersymmetry 

breaking, the gravitino or any standard model superpartner could in principle be the 

LSP. In contrast, in "low-scale" supersymmetry breaking models with M m <C M*, 

such as gauge-mediated supersymmetry breaking models, 

F F 
fn 3/2 = - < m/, m A ~ — — , (12) 
V3M, ; M m 

\/F ~ y/M weak M m <C 10 10 GeV, and the gravitino is necessarily the LSP. 

As with electroweak symmetry breaking, the dynamics of supersymmetry break- 
ing contributes to the energy density of the vacuum, that is, to dark energy. In non- 
supersymmetric theories, the vacuum energy density is presumably naturally A ~ 
instead of its measured value ~ meV 4 , a discrepancy of 10 120 . This is the cosmological 
constant problem. In supersymmetric theories, the vacuum energy density is naturally 
F 2 . For high-scale supersymmetry breaking, one finds A ~ M 2 eak M 2 , reducing the 
discrepancy to 10 90 . Lowering the supersymmetry breaking scale as much as possible 
to F ~ M 2 eak gives A ~ M^ eak , still a factor of 10 60 too big. Supersymmetry there- 
fore eliminates from 1/4 to 1/2 of the fine-tuning in the cosmological constant, a truly 
underwhelming achievement. One must look deeper for insights about dark energy and 
a solution to the cosmological constant problem. 



2.6 Minimal Supergravity 

To obtain detailed information regarding the superpartner spectrum, one must turn to 
specific models. These are motivated by the expectation that the weak- scale supersym- 
metric theory is derived from a more fundamental framework, such as a grand unified 
theory or string theory, at smaller length scales. This more fundamental theory should 
be highly structured for at least two reasons. First, unstructured theories lead to vio- 
lations of low energy constraints, such as bounds on flavor-changing neutral currents 
and CP- violation in the kaon system and in electric dipole moments. Second, the gauge 
coupling constants unify at high energies in supersymmetric theories P and a more fun- 
damental theory should explain this. 

From this viewpoint, the many parameters of weak-scale supersymmetry should be 
derived from a few parameters defined at smaller length scales through renormaliza- 
tion group evolution. Minimal supergravity) 6 1 7 1 8 1 9 1 ^* the canonical model for studies of 
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Fig. 5. Renormalization group evolution of supersymmetric mass parameters. From 
Ref.DU 

supersymmetry phenomenology and cosmology, is defined by 5 parameters: 

m Q , Mi/2, Ai, tan /3, sign (//) , (13) 

where the most important parameters are the universal scalar mass m and the universal 
gaugino mass M 1 / 2 , both defined at the grand unified scale M GUT ~ 2 x 10 16 GeV. In 
fact, there is a sixth free parameter, the gravitino mass 

m 3/2 . (14) 

As noted in Sec. 12.51 the gravitino may naturally be the LSP. It may play an important 
cosmological role, as we will see in Sec.lU For now, however, we follow most of the 
literature and assume the gravitino is heavy and so irrelevant for most discussions. 

The renormalization group evolution of supersymmetry parameters is shown in 
Fig. |3 for a particular point in minimal supergravity parameter space. This figure il- 
lustrates several key features that hold more generally. First, as superpartner masses 
evolve from M G ut to M wea k, gauge couplings increase these parameters, while Yukawa 
couplings decrease them. At the weak scale, colored particles are therefore expected 
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to be heavy, and unlikely to be the LSP. The Bino is typically the lightest gaugino, and 
the right-handed sleptons (more specifically, the right-handed stau tr) are typically the 
lightest scalars. 

Second, the mass parameter is typically driven negative by the large top 
Yukawa coupling. This is a requirement for electroweak symmetry breaking: at tree- 
level, minimization of the electroweak potential at the weak scale requires 

m 2 H -m 2 H t&ri 2 p 1 2 
H = tan^-1 2 mz 

where the last line follows for all but the lowest values of tan (3, which are phenomeno- 
logically disfavored anyway. Clearly, this equation can only be satisfied if m 2 Hu < 0. 
This property of evolving to negative values is unique to rn 2 Hu \ all other mass param- 
eters that are significantly diminished by the top Yukawa coupling also experience a 
compensating enhancement from the strong gauge coupling. This behavior naturally 
explains why SU(2) is broken while the other gauge symmetries are not. It is a beauti- 
ful feature of supersymmetry derived from a simple high energy framework and lends 
credibility to the extrapolation of parameters all the way up to a large mass scale like 
Mgut- 

Given a particular high energy framework, one may then scan parameter space to 
determine what possibilities exist for the LSP. The results for a slice through minimal 
supergravity parameter space are shown in Fig. They are not surprising. The LSP 
is either the the lightest neutralino \ or the right-handed stau tr. In the \ LSP case, 
contours of gaugino-ness 

R x = l a sl 2 + \ a w\ 2 > ( 16 ) 

where 

X = ag{-%B) + a w (-iW) + a^fi A + a &u H u , (17) 

are also shown. The neutralino is nearly pure Bino in much of parameter space, but may 
have a significant Higgsino mixture for m ^ 1 TeV, where Eq. <fl3T> implies ~ Mi. 

There are, of course, many other models besides minimal supergravity. Phenomena 
that do not occur in minimal supergravity may very well occur or even be generic 
in other supersymmetric frameworks. On the other hand, if one looks hard enough, 
minimal supergravity contains a wide variety of dark matter possibilities, and it will 
serve as a useful framework for illustrating many points below. 
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Fig. 6. Regions of the (m , Ml/ 2 ) parameter space in minimal supergravity with A = 
0, tan/3 = 10, and \i > 0. The lower shaded region is excluded by the LEP chargino 
mass limit. The stau is the LSP in the narrow upper shaded region. In the rest of 
parameter space, the LSP is the lightest neutralino, and contours of its gaugino-ness R x 
(in percent) are shown. From Ref. 12 

2.7 Summary 

• Supersymmetry is a new spacetime symmetry that predicts the existence of a new 
boson for every known fermion, and a new fermion for every known boson. 

• The gauge hierarchy problem may be solved by supersymmetry, but requires that 
all superpartners have masses at the weak scale. 

• The introduction of superpartners at the weak scale mediates proton decay at un- 
acceptably large rates unless some symmetry is imposed. An elegant solution, 
i?-parity conservation, implies that the LSP is stable. Electrically neutral super- 
partners, such as the neutralino and gravitino, are therefore promising dark matter 
candidates. 

• The superpartner masses depend on how supersymmetry is broken. In models 
with high-scale supersymmetry breaking, such as supergravity, the gravitino may 
or may not be the LSP; in models with low-scale supersymmetry breaking, the 
gravitino is the LSP. 
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• Among standard model superpartners, the lightest neutralino naturally emerges 
as the dark matter candidate from the simple high energy framework of minimal 
supergravity. 

• Supersymmetry reduces fine tuning in the cosmological constant from 1 part in 
10 120 to 1 part in 10 60 to 10 90 , and so does not provide much insight into the 
problem of dark energy. 

3 Neutralino Cosmology 

Given the motivations described in Sec. |2]for stable neutralino LSPs, it is natural to 
consider the possibility that neutralinos are the dark matterEU^ 1 j n m is section, we 
review the general formalism for calculating thermal relic densities and its implications 
for neutralinos and supersymmetry. We then describe a few of the more promising 
methods for detecting neutralino dark matter. 



Dark matter may be produced in a simple and predictive manner as a thermal relic 
of the Big Bang. The very early universe is a very simple place — all particles are 
in thermal equilibrium. As the universe cools and expands, however, interaction rates 
become too low to maintain this equilibrium, and so particles "freeze out." Unstable 
particles that freeze out disappear from the universe. However, the number of stable 
particles asymptotically approaches a non-vanishing constant, and this, their thermal 
relic density, survives to the present day. 

This process is described quantitatively by the Boltzmann equation 



where n is the number density of the dark matter particle x, H is the Hubble parameter, 
(uav) is the thermally averaged annihilation cross section, and n eq is the x number den- 
sity in thermal equilibrium. On the right-hand side of Eq. (fT8l) . the first term accounts 
for dilution from expansion. The n 2 term arises from processes XX ~^ ff that destroy 
X particles, and the term arises from the reverse process // — > XX* which creates 
X particles. 

It is convenient to change variables from time to temperature, 



3.1 Freeze Out and WIMPs 
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Fig. 7. The co-moving number density Y of a dark matter particle as a function of 
temperature and time. From Ref. |TS| 

where m is the x mass, and to replace the number density by the co-moving number 
density 

n 

n^Y = - , (20) 

s 

where s is the entropy density. The expansion of the universe has no effect on Y, 
because s scales inversely with the volume of the universe when entropy is conserved. 
In terms of these new variables, the Boltzmann equation is 

x dY = n eq (a A v) ( Y 2 \ 
Y cq dx ~ H \Y C \ ) • 

In this form, it is clear that before freeze out, when the annihilation rate is large com- 
pared with the expansion rate, Y tracks its equilibrium value Y cq . After freeze out, Y 
approaches a constant. This constant is determined by the annihilation cross section 
(uav). The larger this cross section, the longer Y follows its exponentially decreasing 
equilibrium value, and the lower the thermal relic density. This behavior is shown in 
Fig. 13 

Let us now consider WIMPs — weakly interacting massive particles with mass and 
annihilation cross section set by the weak scale: m 2 ~ (cr^v ) _1 ~ M^ eak . Freeze out 
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takes place when 

n cq (cr A v) ~ H . (22) 

Neglecting numerical factors, n cq ~ (mT) 3//2 e _m//T for a non-relativistic particle, and 
H ~ T 2 /M*. From these relations, we find that WIMPs freeze out when 

i \ 1/2" 



m 

— ~ In 

T 



(Tfl \ 
TJ 



30 . (23) 



Since \mv 2 = §T, WIMPs freeze out with velocity v ~ 0.3. 

One might think that, since the number density of a particle falls exponentially once 
the temperature drops below its mass, freeze out should occur at T ~ m. This is not 
the case. Because gravity is weak and M* is large, the expansion rate is extremely slow, 
and freeze out occurs much later than one might naively expect. For a m ~ 300 GeV 
particle, freeze out occurs not at T ~ 300 GeV and time t ~ 10~ 12 s, but rather at 
temperature T ~ 10 GeV and time t ~ 10~ 8 s. 

With a little more work,E3 one can find not just the freeze out time, but also the 
freeze out density 

10- 10 GeV ' 
(ctav) 

A typical weak cross section is 



tt x = msY(x = oo) ; ; . (24) 



2 

M ~ ^ 10- 9 GeV- 2 , (25) 

corresponding to a thermal relic density of Vth 2 ~ 0.1. WIMPs therefore naturally 
have thermal relic densities of the observed magnitude. The analysis above has ignored 
many numerical factors, and the thermal relic density may vary by as much as a few 
orders of magnitude. Nevertheless, in conjunction with the other strong motivations for 
new physics at the weak scale, this coincidence is an important hint that the problems 
of electro weak symmetry breaking and dark matter may be intimately related. 



3.2 Thermal Relic Density 

We now want to apply the general formalism above to the specific case of neutralinos. 
This is complicated by the fact that neutralinos may annihilate to many final states: 
//, W + W~, ZZ, Zh, hh, and states including the heavy Higgs bosons H, A, and H ± . 
Many processes contribute to each of these final states, and nearly every supersymmetry 
parameter makes an appearance in at least one process. The full set of annihilation 



17 



diagrams is discussed in Ref. [18j Codes to calculate the relic density are publicly 
available.^ 

Given this complicated picture, it is not surprising that there are a variety of ways to 
achieve the desired relic density for neutralino dark matter. What is surprising, however, 
is that many of these different ways may be found in minimal supergravity, provided 
one looks hard enough. We will therefore consider various regions of minimal super- 
gravity parameter space where qualitatively distinct mechanisms lead to neutralino dark 
matter with the desired thermal relic density. 

3.2.1 Bulk Region 

As evident from Fig. the LSP is a Bino-like neutralino in much of minimal super- 
gravity parameter space. It is useful, therefore, to begin by considering the pure Bino 
limit. In this case, all processes with final state gauge bosons vanish. This follows from 
supersymmetry and the absence of 3-gauge boson vertices involving the hypercharge 
gauge boson. 

The process XX ff through a t-channel sfermion does not vanish in the Bino 
limit. This process is the first shown in Fig. [8] This reaction has an interesting structure. 
Recall that neutralinos are Majorana fermions. If the initial state neutralinos are in an 
5*-wave state, the Pauli exclusion principle implies that the initial state is CP-odd, with 
total spin S = and total angular momentum J = 0. If the neutralinos are gauginos, the 
vertices preserve chirality, and so the final state //has spin S = 1. This is compatible 
with J = only with a mass insertion on the fermion line. This process is therefore 
either P-wave-suppressed (by a factor v 2 ~ 0.1) or chirality suppressed (by a factor 
rrif/M w ). In fact, this conclusion holds also for mixed gaugino-Higgsino neutralinos 
and for all other processes contributing to the // final state. 18 (It also has important 
implications for indirect detection. See Sec. 13.41 ) 

The region of minimal supergravity parameter space with a Bino-like neutralino 
where XX ff yields the right relic density is the (m , M"i/ 2 ) ~ (100 GeV, 200 GeV) 
region shown in Fig. |5J It is called the "bulk region," as, in the past, there was a 
wide range of parameters with m , Mi/ 2 ^ 300 GeV that predicted dark matter within 
the observed range. The dark matter energy density has by now become so tightly 
constrained, however, that the "bulk region" has now been reduced to a thin ribbon of 
acceptable parameter space. 

Moving from the bulk region by increasing m and keeping all other parameters 
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Fig. 8. Three representative neutralino annihilation diagrams. 

fixed, one finds too much dark matter. This behavior is evident in Fig. |9] and not diffi- 
cult to understand: in the bulk region, a large sfermion mass suppresses (cr^f), which 
implies a large fi DM . In fact, sfermion masses not far above current bounds are required 
to offset the P-wave suppression of the annihilation cross section. This is an interest- 
ing fact — cosmology seemingly provides an upper bound on superpartner masses! If 
this were true, one could replace subjective naturalness arguments by the fact that the 
universe cannot be overdosed to provide upper bounds on superpartner masses. 

Unfortunately, this line of reasoning is not airtight even in the constrained frame- 
work of minimal supergravity. The discussion above assumes that \X ~> ff is me only 
annihilation channel. In fact, however, for non-Bino-like neutralinos, there are many 
other contributions. Exactly this possibility is realized in the focus point region, which 
we describe next. 

In passing, it is important to note that the bulk region, although the most straight- 
forward and natural in many respects, is also severely constrained by other data. The 
existence of a light superpartner spectrum in the bulk region implies a light Higgs boson 
mass, and typically significant deviations in low energy observables such as b — > sj 
and (g — 2) M . Current bounds on the Higgs boson mass, as well as concordance be- 
tween experiments and standard model predictions for b — ► 57 and (possibly) (g — 2) M , 
therefore disfavor this region, as can be seen in Fig. |9j For this reason, it is well worth 
considering other possibilities, including the three we now describe. 
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Fig. 9. The bulk and co-annihilation regions of minimal supergravity with A = 0, 
tan /3 = 10 and [i < 0. In the light blue region, the thermal relic density satisfies the 
pre-WMAP constraint 0.1 < Qdm^ 2 < 0.3. In the dark blue region, the neutralino 
density is in the post-WMAP range 0.094 < fi D M^ 2 < 0.129. The bulk region is the 
dark blue region with (m ,M 1/2 ) ~ (100 GeV, 200 GeV). The stau LSP region is 
given in dark red, and the co-annihilation region is the dark blue region along the stau 
LSP border. Current bounds on b — > S7 exclude the green shaded region, and the Higgs 
mass is too low to the left of the nih = 114 GeV contour. From Ref.l20l 

3.2.2 Focus Point Region 

As can be seen in Fig.0 a Bino-like LSP is not a definitive prediction of minimal su- 
pergravity. For large scalar mass parameter m , the Higgsino mass parameter \n\ drops 
to accommodate electroweak symmetry breaking, as required by Eq. (fT5l) . The LSP 
then becomes a gaugino-Higgsino mixture. The region where this happens is called 
the focus point region, a name derived from peculiar properties of the renormaliza- 
tion group equations which suggest that large scalar masses do not necessarily imply 
fine-tuning ™3I 

In the focus point region, the first diagram of Fig. [S] is suppressed by very heavy 
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Fig. 10. Focus point region of minimal supergravity for A = 0, // > 0, and tan (5 as 
indicated. The excluded regions and contours are as in Fig.|6l In the light yellow region, 
the thermal relic density satisfies the pre-WMAP constraint 0.1 < Odm^ 2 < 0.3. In 
the medium red region, the neutralino density is in the post-WMAP range 0.094 < 
Odm^ 2 < 0.129. The focus point region is the cosmologically favored region with 
m > 1 TeV. Updated from Ref.QU 

sfermions. However, the existence of Higgsino components in the LSP implies that 
diagrams like the 2nd of Fig.[U XX ~^ W + W~ through a it-channel chargino, are no 
longer suppressed. This provides a second method by which neutralinos may annihilate 
efficiently enough to produce the desired thermal relic density. The cosmologically 
preferred regions with the right relic densities are shown in Fig. [TOj The right amount 
of dark matter can be achieved with arbitrarily heavy sfermions, and so there is no 
useful cosmological upper bound on superpartner masses, even in the framework of 
minimal supergravity. 

3.2.3 A Funnel Region 

A third possibility realized in minimal supergravity is that the dark matter annihilates to 
fermion pairs through an s-channel pole. The potentially dominant process is through 
the A Higgs boson (the last diagram of Fig. [U), as the A is CP-odd, and so may couple 
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Fig. 11. The A funnel region of minimal supergravity with A = 0, tan/3 = 45, and 
\i < 0. The red region is excluded. The other shaded regions have f2 DM ^ 2 < 0.1 
(yellow), 0.1 < DM /i 2 < 0.3 (green), and 0.3 < Q DM h 2 < 1 (blue). FromRef.EU 

to an initial S'-wave state. This process is efficient when 2m x « m^. In fact, the 
A resonance may be broad, extending the region of parameter space over which this 
process is important. 

The A resonance region occurs in minimal supergravity for tan (3 ^ 4J 24 I 25 I and is 
shown in Fig. [TT] Note that the resonance is so efficient that the relic density may be 
reduced too much. The desired relic density is therefore obtained when the process is 
near resonance, but not exactly on it. 

3.2.4 Co-annihilation Region 

Finally, the desired neutralino relic density may be obtained even if \X annihilation 
is inefficient if there are other particles present in significant numbers when the LSP 
freezes out. The neutralino density may then be brought down through co-annihilation 
with the other speciesP^ZI Naively, the presence of other particles requires that they 
be mass degenerate with the neutralino to within the temperature at freeze out, T « 
m x /30. In fact, co-annihilation may be so enhanced relative to the P-wave-suppressed 
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XX annihilation cross section that co-annihilation may be important even with mass 
splittings much larger than T. 

The co-annihilation possibility is realized in minimal supergravity along the f LSP 
- x LSP border. (See Fig. 0) Processes such as %f — ► t* — ► tj are not P-wave 
suppressed, and so enhance the XX annihilation rate substantially. There is therefore 
a narrow finger extending up to masses m x ~ 600 GeV with acceptable neutralino 
thermal relic densities. 



$ x = 6.6 x 10 4 cm- s-i —g— - - - . (26) 



3.3 Direct Detection 

If dark matter is composed of neutralinos, it may be detected directly, that is, by look- 
ing for signals associated with its scattering in ordinary matter. Dark matter veloc- 
ity and spatial distributions are rather poorly known and are an important source of 
uncertainty ! 28 | 29 | 30 |3 - jJ3 ^ 1 A common assumption is that dark matter has a local energy 
density of p x = 0.3 GeV/cm 3 with a velocity distribution characterized by a velocity 
v ~ 220 km/s. Normalizing to these values, the neutralino flux is 

,4 _-a „-i P X 100 GeV 

0.3GeV/cm 3 m x 220 km/s ' 

Such values therefore predict a substantial flux of halo neutralinos in detectors here on 
Earth. 

The maximal recoil energy from a WIMP scattering off a nucleus iV is 

2rrii.m N „ 

= 7 v 2 ~ 100 keV . (27) 

[m x + m N ) 

With such low energies, elastic scattering is the most promising signal at present, al- 
though the possibility of detecting inelastic scattering has also been discussed. As we 
will see below, event rates predicted by supersymmetry are at most a few per kilo- 
gram per day. Neutralino dark matter therefore poses a serious experimental challenge, 
requiring detectors sensitive to extremely rare events with low recoil energies. 

Neutralino-nucleus interactions take place at the parton level through neutralino- 
quark interactions, such as those in Fig. [12| Because neutralinos now have velocities 
v ~ 10~ 3 , we may take the non-relativistic limit for these scattering amplitudes. In this 
limit, only two types of neutralino-quark couplings are non- vanishing. 33 The interac- 
tion Lagrangian may be parameterized as 

£ = E (ufxi^xqirfq + afxxqq) ■ (28) 

q=u,d,s,c,b,t 
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Fig. 12. Feynman diagrams contributing to XQ — > XQ scattering. 

The first term is the spin-dependent coupling, as it reduces to S x ■ Sn in the non- 
relativistic limit. The second is the spin-independent coupling. All of the supersymme- 
try model dependence is contained in the parameters and a^ 1 . The t-channel Higgs 
exchange diagram of Fig, [recontributes solely to a^ 1 , while the s-channel squark dia- 
gram contributes to both a q D and a^ 1 . 

For neutralinos scattering off protons, the spin-dependent coupling is dominant. 
However, the spin-independent coupling is coherent and so greatly enhanced for 
heavy nuclei, a fact successfully exploited by current experiments. As a result, spin- 
independent direct detection is currently the most promising approach for neutralino 
dark matter, and we focus on this below. 

Given the parameters a^ 1 , the spin-independent cross section for xN scattering is 



2 SI 2 

71 q 



Z^f Tq + (A-Z)^fl 
m q q m q q 



(29) 



where 



m x + m N 

is the reduced mass of the system, Z and A are the atomic number and weight of 
the nucleus, respectively, and 

fP ,n _ (p,n\m q qq\p,n} 

JT q — P-U 

are constants quantifying what fraction of the nucleon's mass is carried by quark q. For 
the light quarks J*^ 

ft = 0.020 ±0.004 ft =0.014 ±0.003 
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/£ = 0.026 ± 0.005 f% 4 = 0.036 ± 0.008 

/£ = 0.118 ±0.062 ft = 0.118 ± 0.062 . (32) 



The contribution from neutralino-gluon couplings mediated by heavy quark loops may 
be included by taking f"^ = = £(1 - - f% - BB 

The number of dark matter scattering events is 

TV = N N T^a N v (33) 
m x 

3.4 x iq-6 M ^ T Px 100 GeV I iM °JL , (34) 

kg day 0.3 GeV/cm 3 m x 220km/s m 2 p 10~ 6 pb 

where N N is the number of target nuclei, T is the experiment's running time, M D is the 
mass of the detector, and the proton scattering cross section u v has been normalized to a 
near-maximal supersymmetric value. This is a discouragingly low event rate. However, 
for a detector with a fixed mass, this rate is proportional to n 2 N A. For heavy nuclei with 
A ~ m x /m p , the event rate is enhanced by a factor of ~ A 3 , providing the strong 
enhancement noted above. 

Comparisons between theory and experiment are typically made by converting all 
results to proton scattering cross sections. In Fig. [T31 minimal supergravity predic- 
tions for spin-independent cross sections are given. These vary by several orders of 
magnitude. In the stau co-annihilation region, these cross sections can be small, as the 
neutralino is Bino-like, suppressing the Higgs diagram, and squarks can be quite heavy, 
suppressing the squark diagram. However, in the focus point region, the neutralino is 
a gaugino-Higgsino mixture, and the Higgs diagram is large. Current and projected 
experimental sensitivities are also shown in Fig. [T3j Current experiments are just now 
probing the interesting parameter region for supersymmetry, but future searches will 
provide stringent tests of some of the more promising minimal supergravity predic- 
tions. 

The DAMA collaboration has reported evidence for direct detection of dark matter 
from annual modulation in scattering rates. 37 The favored dark matter mass and proton 
spin-independent cross section are shown in Fig.fTH By comparing Figs. [T3l and [T?l 
one sees that the interaction strength favored by DAMA is very large relative to typical 
predictions in minimal supergravity. Such cross sections may be realized in less re- 
strictive supersymmetry scenarios. However, more problematic from the point of view 
of providing a supersymmetric interpretation is that the experiments EDELWEISS^ 
and CDMS 39 have also searched for dark matter with similar sensitivities and have not 
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Fig. 13. Spin-independent neutralino-proton cross sections for minimal supergravity 
models with A = 0, fi > 0. The colors correspond to various values of tan/3 
in the range 10 < tan/3 < 55. Points with the preferred thermal relic density 
0.094 < f^DM^ 2 < 0.129 are highlighted with enlarged circles, and those in the focus 
point and co-annihilation regions are indicated. Estimated reaches of current (CDMS, 
EDELWEISS, ZEPLIN1, DAMA), near future (CDMS2, EDELWEISS2, ZEPLIN2, 
CRESST2), and future detectors (GENIUS, ZEPLIN4,CRYOARRAY) are given by the 
solid, dark dashed, and light dashed contours, respectively. From Ref.l36l 

found signals. Their exclusion bounds are also given in Fig.[T3|* Given standard halo 
and neutralino interaction assumptions, these data are inconsistent at a high level. 

Non-standard halo models and velocity distribution^ 1 ^ 1 and non-standard and gen- 
eralized dark matter interaction j 43 | 44 | 45 | ^l have been considered as means to bring con- 
sistency to the experimental picture. The results are mixed. Given the current status 
of direct detection experiments, a supersymmetric interpretation is at best premature. 
It is worth noting, however, that the current results bode well for the future, as many 
well-motivated supersymmetry models predict cross sections not far from current sen- 
sitivities. 

Recent data from CDMS in the Soudan mine has pushed the discrepancy to even greater levelsP 
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WIMP Mass [GeV/c 2 ] 

Fig. 14. Regions of dark matter mass and spin-independent proton scattering cross 
sections. The shaded region is the 3<r favored region from DAMA. The dot-dashed 
line is the exclusion contour from EDELWEISS, and the thick solid black line is the 
exclusion contour from CDMS. From Ref. 39. 

3.4 Indirect Detection 

After freeze out, dark matter pair annihilation becomes greatly suppressed. However, 
after the creation of structure in the universe, dark matter annihilation in overdense re- 
gions of the universe may again become significant. Dark matter may therefore be de- 
tected indirectly: pairs of dark matter particles annihilate somewhere, producing some- 
thing, which is detected somehow. There are a large number of possibilities. Below we 
briefly discuss three of the more promising signals. 

3.4.1 Positrons 

Dark matter in our galactic halo may annihilate to positrons, which may be detected 
in space-based or balloon-borne experiments! 47 1 48 1 49 1 50 1 51 1 52 ' (Anti-proton s' 53 1 54 1 5 - 5 ^ and 
anti-deuteriuirPZIhave also been suggested as promising signals.) 

The positron background is most likely to be composed of secondaries produced 
in the interactions of cosmic ray nuclei with interstellar gas, and is expected to fall as 
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~ E~+ s . At energies below 10 GeV, there are also large uncertainties in the back- 
ground.^^ The most promising signal is therefore hard positrons from XX annihila- 
tion. 

Unfortunately, the monoenergetic signal XX ~> e + e _ is extremely suppressed. As 
noted above, XX ~ > // is either P-wave suppressed or chirality suppressed. At present 
times, as opposed to during freeze out, P-wave suppression is especially severe, since 
v 2 ~ 1CT 6 , and so direct annihilation to positrons is effectively absent. ^ The positron 
signal therefore results from processes such as XX ~ * W + W~ followed by W + — > e + u, 
and is a continuum, not a line, at the source. 

To obtain the positron energy distribution we would observe, the source energy dis- 
tribution must be propagated through the halo to us. The resulting differential positron 

fiux ysa 2 

^= = ^\J2°i vB U I dE fi(E ) G(E Q ,E) , (35) 
aildh mi J 

x l 

where p x is the local neutralino mass density, the sum is over all annihilation channels, 
and B l e+ is the branching fraction to positrons in channel i. The source function f(E ) 
gives the initial positron energy distribution from neutralino annihilation. G(E , E) is 
the Green's function describing positron propagation in the galaxy^ and contains all 
the halo model dependence. 

Three sample positron spectra are given in Fig. [15] For all of them, E 2 d$/dE 
peaks at energies E ~ m x /2. These signals are all well below background. However, 
a smooth halo distribution has been assumed. For clumpy halos, which are well within 
the realm of possibility, the signal may be enhanced significantly. In the next few years, 
both PAMELA, a satellite detector, and AMS-02, an experiment to be placed on the 
International Space Station, will provide precision probes of the positron spectrum. 
These experiments and other recently completed experiments are listed in Table [lj 

Finally, the High Energy Antimatter Telescope (HEAT) experiment, a balloon- 
borne magnetic spectrometer, has found evidence for an excess of positrons at energy 
~ 8 GeV in data from 1994/95™ 1 and 2000P The observed bump in the positron 
fraction N e + / (N e + + N e - ) is not naturally obtained by neutralino dark matter for two 
reasons. First, as noted above, for a smooth halo, the annihilation cross sections that 
produce the desired relic density predict positron fluxes that are far too low to explain 
the observed excess. In principle, this objection may be overcome by a sufficiently 

tNote that this suppression is rather special, in that it follows from the Majorana nature of neutralinos; it 
is absent for other types of dark matter, such as dark matter with spin 1 .HUSO 
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Fig. 15. The differential positron flux for three minimal supergravity models. The 
curves labeled C and HEMN are background models from Ref. |52| From Ref. 



Table 1. Recent and planned e + detector experiments. We list each experiment's start 
date, duration, geometrical acceptance in cm 2 sr, maximal E e + sensitivity in GeV, and 
(expected) total number of e + detected per GeV at E e + = 100 GeV. From Ref. 95 



Experiment 



Type 



Date Duration Accept. 



j^max dN 



dE 



(100) 



HEAT94/95 
CAPRICE94/98 
PAMELA 
AMS-02 



Balloon 
Balloon 
Satellite 
ISS 



1994/95 
1994/98 



29/26 hr 
18/21 hr 
3 yr 
3 yr 



495 
163 
20 
6500 



50 
10/30 
200 
1000 



0.7 
250 



clumpy halo. Second, neutralino annihilation produces positrons only through cas- 
cades, resulting in a smooth positron energy distribution. This is an inevitable con- 
sequence of the Majorana nature of neutralinos. Nevertheless, even the addition of a 
smooth component from neutralino annihilation may improve the fit to data, and the 
possibility of a supersymmetric explanation for the HEAT anomaly has been explored 
in a number of studies E2I85ESI 

Two "best fit" results from Ref.l65lare shown in Fig.[l6j In this study, the ignorance 
of subhalo structure is parameterized by a constant B s , an overall normalization factor 
that enhances the positron flux relative to what would be expected for a smooth halo. 
As can be seen in Fig. [T^J both spectra give improved fits to the data. They require 
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Fig. 16. Two positron spectra for which contributions from neutralino annihilation im- 
prove the fit to HEAT data. In each case, the contribution from neutralino annihilation 
has been enhanced by a factor B s ~ 100 relative to the prediction for a smooth halo. 
From Ref .1651 

substantial boost factors, however, with B s ~ 100. Such large boost factors may be 
disfavored by models of halo formation.^ 

3.4.2 Photons 

Dark matter in the galactic center may annihilate to photons, which can be de- 
tected in atmospheric Cherenkov telescopes on the ground, or in space-based detec- 
tors .EESESEaznZlZSEll (Photons from the galactic haloF^or even from extra-galactic 
source^ have also been considered.) 

The main source of photons is from cascade decays of other primary annihilation 
products. A line source from loop-mediated processes such as \X ~^ 70 * 78 | 79 | 8U * and 
XX ~> 7^^is possible,^ but is typically highly suppressed.^ 

The differential photon flux along a direction that forms an angle ip with respect to 
the direction of the galactic center is 

5MB = ?i ff<0 ei/ dl ' (36) 

* A. 

where the sum is over all annihilation channels, p is the neutralino mass density, and 
the integral is along the line of sight. All of the halo model dependence is isolated in 
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Fig. 17. Integral photon fluxes $ 7 (£ , t h r ) as a function of threshold energy 
i?thr for A = 0, ji > 0, m t = 174 GeV, and halo parameter J = 
500. The four models have relic density VL x h 2 « 0.15, and are specified by 
(tan (3, mo, M 1/2 , m x , R x ) = (10, 100, 170, 61, 0.93) (dotted), (10, 1600, 270, 97, 0.77) 
(dashed), (10, 2100, 500, 202, 0.88) (dot-dashed), and (50, 1000, 300, 120, 0.96) (solid), 
where all masses are in GeV. Point source flux sensitivity estimates for several gamma 
ray detectors are also shown. From Ref. 



the integral. Depending on the dumpiness or cuspiness of the halo density profile, this 
integral may vary by as much as 5 orders of magnitude.^ 

The integrated photon signal for 4 representative minimal supergravity models is 
given in Fig. [l7| A moderate halo profile is assumed. Experiments sensitive to such 
photon fluxes are listed in Table|2j and their sensitivities are given in Fig.fTTl 



3.4.3 Neutrinos 

When neutralinos pass through astrophysical objects, they may be slowed below es- 
cape velocity by elastic scattering. Once captured, they then settle to the center, where 
their densities and annihilation rates are greatly enhanced. While most of their annihi- 
lation products are immediately absorbed, neutrinos are not. Neutralinos may there- 
fore annihilate to high energy neutrinos in the cores of the EartrP™OTH and 
g UI j85 | 87 | 88 | 89 | 90 | 9i | 92 | 93 | 94 | an( ^ ^ e detected on Earth in neutrino telescopes. 

The formalism for calculating neutrino fluxes from dark matter annihilation is com- 
plicated but well developed. (See Ref. 16 for a review.) In contrast to the previous 
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Table 2. Some of the current and planned 7 ray detector experiments with sensitivity 
to photon energies 10 GeV < E 1 1z 300 GeV. We list each experiment's start date and 
expected E 1 coverage in GeV. The energy ranges are approximate. For experiments 
constructed in stages, the listed threshold energies will not be realized initially. From 
Ref.El 



Experiment 


Type 


Date 


£L Range 


EGRET 


Satellite 


l 991-2000 


0.02-30 


STACEE 


ACT array 


1998 


20-300 


CELESTE 


ACT array 


1998 


20-300 


ARGO-YBJ 


Air shower 


2001 


100-2,000 


MAGIC 


ACT 




10-1000 


AGILE 


Satellite 




0.03-50 


HESS 


ACT array 




10-1000 


AMS/7 


Space station 




0.3-100 


CANGAROO III 


ACT array 




30-50,000 


VERITAS 


ACT array 




50-50,000 


GLAST 


Satellite 




0.1-300 



two indirect detection examples, neutrino rates depend not only on annihilation cross 
sections, but also on xN scattering, which determines the neutralino capture rate in the 
Sun and Earth. 

As with positrons, XX ~ > vv is helicity-suppressed, and so neutrinos are produced 
only in the decays of primary annihilation products. Typical neutrino energies are then 
E v ~ |m x to \m x , with the most energetic spectra resulting from WW, ZZ, and, 
to a lesser extent, rf. After propagating to the Earth's surface, neutrinos are detected 
through their charged-current interactions. The most promising signal is from upward- 
going muon neutrinos that convert to muons in the surrounding rock, water, or ice, 
producing through-going muons in detectors. The detection rate for such neutrinos is 
greatly enhanced for high energy neutrinos, as both the charged-current cross section 
and the muon range are proportional to E v . 

The most promising source of neutrinos is the core of the Sun. Muon flux rates 
from the Sun are presented in Fig.[T8j Fluxes as large as 1000 km -2 s -1 are possible. 
Past, present, and future neutrino telescopes and their properties are listed in Table |3] 
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Fig. 18. Muon flux from the Sun in km 2 yr 1 for v = 270 km/s and p x 
0.3 GeV/cm 3 . FromRef.| 



Comparing Fig. [T8] with Table we find that present limits do not significantly con- 
strain the minimal supergravity parameter space. However, given that the effective area 
of neutrino telescope experiments is expected to increase by 10 to 100 in the next few 
years, muon fluxes of order 10-100 km~ 2 yr -1 may be within reach. 



3.5 Summary 

Neutralinos are excellent dark matter candidates. The lightest neutralino emerges nat- 
urally as the lightest supersymmetric particle and is stable in simple supersymmetric 
models. In addition, the neutralino is non-baryonic, cold, and weakly-interacting, and 
so has all the right properties to be dark matter, and its thermal relic density is naturally 
in the desired range. 

Current bounds on Odm are already highly constraining. Although these constraints 
do not provide useful upper bounds on supersymmetric particle masses, they do restrict 
supersymmetric parameter space. In minimal supergravity, the cosmologically pre- 
ferred regions of parameter space include the bulk, focus point, A funnel, and stau 
coannihilation regions. 

Neutralinos may be detected either directly through their interactions with ordinary 
matter or indirectly through their annihilation decay products. Null results from direct 
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Table 3. Current and planned neutrino experiments. We list also each experiment's 
start date, physical dimensions (or approximate effective area), muon threshold energy 
_E* hr in GeV, and 90% CL flux limits for the Sun $® in km" 2 yr" 1 for half-cone angle 
6 ps 15° when available. From Ref.l95l 



Experiment 


Type 


Date 


Dimensions 


£thr 

A* 


$© 


Baksan 


Ground 


1978 


17 x 17 x 11 m 3 


1 


7.6 x 10 3 


Kamiokande 


Ground 


1 QCJ 


~ lou m 


1 

3 


1 I X 1U 


MACRO 


Ground 


1989 


12 x 77 x 9 m 3 


2 


6.5 x 10 3 


S uper- Kamiokande 


Ground 


1996 


~ 1200 m 2 


1.6 


5.0 x 10 3 


Baikal NT-96 


Water 


1996 


~ 1000 m 2 


10 




AMANDA B- 10 


Under-ice 


1997 


~ 1000 m 2t 


~ 25 




Baikal NT-200 


Water 


1998 


~ 2000 m 2 






AMANDA II 


Ice 


2000 


~ 3 x 10 4 m 2 


~ 50 




NESTOR § 


Water 




~ 10 4 m 2 * 


few 




ANTARES 


Water 




~2x 10 4 m 2 * 


~ 5-10 




IceCube 


Ice 




~ 10 6 m 2 







t Hard spectrum, m x = 100 GeV. § One tower. 1 ~ 100 GeV. 

and indirect dark matter searches are not yet very constraining. Future sensitivities of 
various particle physics and dark matter detection experiments are shown in Fig. [191 
The sensitivities assumed, and experiments likely to achieve these sensitivities in the 
near future, are listed in Table |4j 

Several interesting features are apparent. First, traditional particle physics and dark 
matter searches, particularly indirect detection experiments, are highly complementary. 
Second, at least one dark matter experiment is predicted to see a signal in almost all 
of the cosmologically preferred region. This illustration is in the context of minimal 
supergravity, but can be expected to hold more generally. The prospects for neutralino 
dark matter discovery are therefore promising. 

4 Gravitino Cosmology 

In Sec. |3j we largely ignored the gravitino. In this Section, we will rectify this omis- 
sion. Although gravitino interactions are highly suppressed, gravitinos may have im- 
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m (GeV) 

Fig. 19. Estimated reaches of various high-energy collider and low-energy precision 
searches (black), direct dark matter searches (red), and indirect dark matter searches 
(blue) in the next few years for minimal supergravity with A = 0, tan/9 = 10, and 
\i > 0. The excluded green regions are as in Fig.|6j The blue (yellow) shaded region 
has 0.1 < Qdm^ 2 < 0.3 (0.025 < Qdm^ 2 < 1)- The regions probed extend the curves 
toward the excluded green regions. From Ref.l95l 

plications for many aspects of cosmology, including Big Bang nucleosynthesis (BBN), 
the cosmic microwave background, inflation, and reheating. Gravitino cosmology is in 
many ways complementary to neutralino cosmology, providing another rich arena for 
connections between microscopic physics and cosmology. 

4.1 Gravitino Properties 

The properties of gravitinos may be systematically derived by supersymmetrizing the 
standard model coupled to gravity. Here we will be content with highlighting the main 
results. 

In an exactly supersymmetric theory, the gravitino is a massless spin 3/2 particle 
with two degrees of freedom. Once supersymmetry is broken, the gravitino eats a spin 
1/2 fermion, the Goldstino of supersymmetry breaking, and becomes a massive spin 
3/2 particle with four degrees of freedom. As noted in Sec. 12.51 the resulting gravitino 
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Table 4. Constraints on supersymmetric models used in Fig. \T9\ We also list experi- 
ments likely to reach these sensitivities in the near future. From Ref. |55J 
Observable Type Bound Experiment(s) 



X ± X° 

Muon MDM 

""proton 

v from Earth 
v from Sun 
7 (gal. center) 
7 (gal. center) 



Collider 
Low energy 
Low energy 
Direct DM 
Indirect DM 
Indirect DM 
Indirect DM 



See Refs. l9"6T9"7l9"%l 

\AB(B -»■ X s j)\ < 1.2 x 10~ 4 
« SY | < 8 x 10- 10 

i fi i 

Fig. E] 

yr 



$® < 100 km" 2 — - 1 



$° < 100 km~ 2 yr" 1 



$J1) < 1.5 x io- 



10 cm 2 s 1 



Indirect DM $ 7 (50) < 7 x 10~ 12 cm~ 2 s" 1 



e + cosmic rays Indirect DM (S/B) majX < 0.01 



Tevatron Run II 
BaBar, BELLE 
Brookhaven E821 
CDMS2, CRESST2 
AMANDA 
AMANDA 
GLAST 
HESS, MAGIC 
AMS-02 



mass is 

rn - = 



m G = 7^77 > ( 37 ) 



where M, = (SttG^)- 1 / 2 ~ 2.4 x 10 18 GeV is the reduced Planck mass. Gravitinos 
couple standard model particles to their superpartners through gravitino-gaugino-gauge 
boson interactions 

L = - si:^ [r ' 7l7 ^ p ' (38) 

and gravitino-sfermion-fermion interactions 

j dJfrf-fG*. (39) 



In models with high-scale supersymmetry breaking, such as conventional super- 
gravity theories, F ~ M wcak M t , as explained in Sec. 12.51 The gravitino mass is there- 
fore of the order of the other superpartner masses, and we expect them all to be in the 
range ~ 100 GeV— ITeV. The gravitino's effective couplings are ~ E/M*, where E is 
the energy of the process. The gravitino's interactions are therefore typically extremely 
weak, as they are suppressed by the Planck scale. 

We will focus on theories with high-scale supersymmetry breaking in the follow- 
ing discussion. Note, however, that in theories with low-scale supersymmetry break- 
ing, the gravitino may be much lighter, for example, as light as ~ eV in some simple 
gauge-mediated supersymmetry breaking models. The gravitino's interactions through 
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its Goldstino components may also be much stronger, suppressed by F/M weak rather 
than M*. For a summary of gravitino cosmology in such scenarios, see Ref.l99l 



4.2 Thermal Relic Density 

If gravitinos are to play a cosmological role, we must first identify their production 
mechanism. There are a number of possibilities. Given our discussion of the neu- 
tralino thermal relic density in Sec. a natural starting place is to consider gravitino 
production as a result of freeze out from thermal equilibrium. At present, the grav- 
itino coupling E/M* is a huge suppression. However, if we extrapolate back to very 
early times with temperatures T ~ M*, even gravitational couplings were strong, and 
gravitinos were in thermal equilibrium, with = n eq . Once the temperature drops 
below the Planck scale, however, gravitinos quickly decouple with the number density 
appropriate for relativistic particles. Following decoupling, their number density then 
satisfies Uq oc R~ 3 oc T 3 . This has the same scaling behavior as the background photon 
number density, however, and so we expect roughly similar number densities now. 

If such gravitinos are stable, they could be dark matter. In fact, the first supersym- 
metric dark matter candidate proposed was the gravitinoP^ However, the overclosure 
bound implies 

tt d < 1 =>■ m 6 < 1 keV . (40) 

This is not surprising — relic neutrinos have a similar density, and the overclosure 
bound on their mass is similar. 

On the other hand, gravitinos may be unstable!^ This may be because imparity 
is broken, or because the gravitino is not the LSP. In this case, there is no bound from 
overclosure, but there are still constraints. In particular, the gravitino's decay products 
will destroy the successful predictions of BBN for light element abundances if the decay 
takes place after BBN. In the case where decay to a lighter supersymmetric particle is 
possible, we can estimate the gravitino lifetime to be 



Ml 

— =- ~ 0.1 yr 

m% 



100 GeV 8 



m G 



(41) 



Requiring gravitino decays to be completed before BBN at t ~ 1 s implies^' 

m d > 10 TeV . (42) 



In both cases, the required masses are incompatible with the most natural expecta- 
tions of conventional supergravity theories. Gravitinos may, however, be a significant 
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component of dark matter if they are stable with mass ~ keV. Such masses are pos- 
sible in low-scale supersymmetry breaking scenarios, given an appropriately chosen 
supersymmetry-breaking scale F. 



4.3 Production during Reheating 

In the context of inflation, the gravitino production scenario of Sec. 14.21 is rather un- 
natural. Between the time when T ~ M* and now, we expect the universe to inflate, 
which would dilute any gravitino relic thermal density. Inflation does provide another 
source for gravitinos, however. Specifically, following inflation, we expect an era of 
reheating, during which the energy of the inflaton potential is transferred to standard 
model particles and superpartners, creating a hot thermal bath in which gravitinos may 



After reheating, the universe is characterized by three hierarchically separated rates: 
the interaction rate of standard model particles and their superpartners with each other, 
ctsm^; the expansion rate, H; and the rate of interactions involving one gravitino, a^n. 
Here n is the number density of standard model particles. After reheating, the universe 
is expected to have a temperature well below the Planck scale, but still well above 
standard model masses. These rates may then be estimated by dimensional analysis, 



The picture that emerges, then, is that after reheating, there is a thermal bath of 
standard model particles and their superpartners. Occasionally these interact to pro- 
duce a gravitino through interactions like gg — > gG. The produced gravitinos then 
propagate through the universe essentially without interacting. If they are stable, as we 
will assume throughout this section, they contribute to the present dark matter density. 

To determine the gravitino abundance, we turn once again to the Boltzmann equa- 
tion 



In this case, the source term nL arises from interactions such as gg — > gG. In contrast 
to our previous application of the Boltzmann equation in Sec. 13.11 however, here the 
n 2 sink term, originating from interactions such as gG — > gg, is negligible. Changing 
variables as before with t — > T and n — > Y = n/s, we find 



beproducedM"^ 1 " 4 ""^ 1 '^ 



and we find 




(43) 




(44) 



dY 
~dT 




(45) 
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process i 


\KA 12 / 9 2 _L m i ^ 
\ M i\ lM^[} + ^k ; ) 


A 


g a + g b -> 5 C + G 


A{s + 2t + 2^)\f abc \ 2 


B 


g a + g b ^g c + G 


-A(t + 2s + 2 s ^)\f abc \ 2 


C 


Qi + g a ^ gj + G 


2s\T^ 


D 


g a + Qi^ qj + G 


-2t\T^\ 2 


E 


9t + q 3 ^g a + G 


-2t\T?f 


F 


~g a + ~g b ^g c + G 


q (s +st+t J rabc\2 
° st(s+t) \J I 


G 


ft + g a ~^q 3 + G 


j / \ lm/i lO 


H 


ft + 9 a -^qj + G 


-2(t + 2s + 2^)|T/.| 2 


I 


ft + qj^g a + G 


"4(t + f)|^| 2 


J 


ft + ?j : -» g a + G 


2(s + 2t + 2f)|T? j | 2 



Fig. 20. Processes contributing to gravitino production after reheating. From Ref. 11071 

The right-hand side is independent of T, since n oc T 3 , H oc T 2 and s oc T 3 . We 
thus find an extremely simple relation — the gravitino relic number density is linearly 
proportional to the reheat temperature T R . 

The constant of proportionality is the gravitino production cross section. The lead- 
ing 2 — > 2 QCD interactions have been included in Ref. 11071 These are listed in 
Fig. El With these results, the gravitino relic density can be determined as a function 
of reheating temperature Tr and gravitino mass. The results are given in Fig.|2T] For 
gravitino mass vtlq ~ 100 GeV, the constraint on Odm requires reheating temperature 
T R ^ 10 10 GeV, providing a bound on the inflaton potential. Of course, if this bound 
is nearly saturated, gravitinos produced after reheating may be a significant component 
of dark matter. 

4.4 Production from Late Decays 

A third mechanism for gravitino production is through the cascade decays of other 
supersymmetric particles. If the gravitino is not the LSR cascade decays will bypass 
the gravitino, given its highly suppressed couplings. However, as discussed in Sec. 12.51 
the gravitino may be the LSR even in high-scale supersymmetry breaking models. If 
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Fig. 21. The gravitino relic abundance VL G h 2 as a function of reheating temperature T R 
for various gravitino masses and gluino mass m~ g = 700 GeV. From Ref. 11071 



the gravitino is the LSP, all cascades will ultimately end in a gravitino. 

An alternative gravitino dark matter scenario is therefore the following! 1U8 | 1U9 1 As- 
sume that the gravitino is the LSP and stable. To separate this scenario from the 
previous two, assume that inflation dilutes the primordial gravitino density and the 
universe reheats to a temperature low enough that gravitino production is negligible. 
Because the gravitino couples only gravitationally with all interactions suppressed by 
the Planck scale, it plays no role in the thermodynamics of the early universe. The 
next-to-lightest supersymmetric particle (NLSP) therefore freezes out as usual; if it is 
weakly-interacting, its relic density will be near the desired value. However, much later, 
after 

M 2 

~ 10 5 s - 10 8 s , (46) 



T ~ 



^weak 



the NLSP decays to the gravitino LSP. The gravitino therefore becomes dark matter 
with relic density 



n 



G = ^NLSP • (47) 

T^NLSP 

The gravitino and NLSP masses are naturally of the same order in theories with high- 
scale supersymmetry breaking. Gravitino LSPs may therefore form a significant relic 
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component of our universe, inheriting the desired relic density from WIMP decay. 
In contrast to the previous two production mechanisms, the desired relic density is 
achieved naturally without the introduction of new energy scales. 

Given our discussion in Sec. 14.21 the decay time of Eq. d4"6"l ). well after BBN, should 
be of concern. In the present case, the decaying particle is a WIMP and so has a density 
far below that of a relativistic particle. (Recall Fig. |7J) However, one must check to 
see if the light element abundances are greatly perturbed. In fact, for some weak-scale 
NLSP and gravitino masses they are, and for some they aren'tP^^ We discuss this 
below, along with other constraints on this scenario. 

Models with weak-scale extra dimensions also provide a similar dark matter par- 
ticle in the form of Kaluza-Klein gravitons) 108 ' 114 1 with Kaluza-Klein gauge bosons or 
leptons playing the role of the decaying WEVIPP^ Because such dark matter can- 
didates naturally preserve the WIMP relic abundance, but have interactions that are 
weaker than weak, they have been named superweakly-interacting massive particles, or 
"superWIMPs."* 

We see now that our discussion in Sec. |3] of WIMP dark matter was only valid for 
the "half" of parameter space where m 3 / 2 > mLSP- When the gravitino is the LSP, there 
are number of new implications of supersymmetry for cosmology. For example, the "f 
LSP" region is no longer excluded by searches for charged dark matterp ^ |109 | 115 | lx ^ 
as the f is no longer stable, but only metastable. There is therefore the possibility of 
stable heavy charged particles appearing in collider detectors) 117 1 118 1 Further, regions 
with too much dark matter are no longer excluded, because the gravitino dark matter 
density is reduced by m^/mNLSP relative to the NLSP density. As we will discuss be- 
low, the late decays producing gravitinos may have detectable consequences for BBN 
and the cosmic microwave background. Astrophysical signatures in the diffuse pho- 
ton spectrum, 108 the ionization of the universep^ and the suppression of small scale 
structure 1 ^ are also interesting possibilities. 

4.5 Detection 

If gravitinos are the dark matter, all direct and indirect searches for dark matter are 
hopeless, because all interaction cross sections and annihilation rates are suppressed by 
the Planck scale. Instead, one must turn to finding evidence for gravitino production in 
the early universe. In the case of gravitinos produced at T ~ M* or during reheating, 

different dark matter candidate that also predicts late decays 
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the relevant physics is at such high energy scales that signals are absent, or at least 
require strong theoretical assumptions. In the case of production by late decays, how- 
ever, there are several possible early universe signals. We consider a few of these in 
this Section. 



4.5.1 Energy Release 

If gravitinos are produced by late decays, the relevant reaction is NLSP — > G + S, 
where S denotes one or more standard model particles. Because the gravitino is 
essentially invisible, the observable consequences rely on finding signals of S pro- 
duction in the early universe. Signals from late decays have been considered in 
Refs. I121I122I123I124I125I126I127I128I129I130I1311 In principle, the strength of 
these signals depends on what S is and its initial energy distribution. It turns out, how- 
ever, that most signals depend only on the time of energy release, that is, the NLSP's 
lifetime r, and the average total electromagnetic or hadronic energy released in NLSP 
decay. 

Here we will consider two possible NLSPs: the photino and the stau. In the photino 

case, 



r(7 - iG) 



1 



48ttM2 m% 

* G 



m 
m 



2 
G 
2 

7 J 



1 + 3-g 



(48) 



In the limit Am 



rrij — uiq <C m^, the decay lifetime is 

^ 7 rlOOGeV 13 

r(7^ 7 G) « 1.8 x 10 7 s 



Am 

independent of the overall superpartner mass scale. For the stau case, 

,2l4 



(49) 



tG) 



1 



48ttM ^ m% 

* G L 



i 



(50) 



In the limit Am 



rrif 



T(T 



^ m G' me decay lifetime is 

s 100 GeV^ 4 
tG) w 3.6 x 10 s s 



Am 



rn 



G 



1 TeV 



(51) 



The electromagnetic energy release is conveniently written in terms of 



,EM 



NLSP 



(52) 



where e E M is the initial electromagnetic energy released in each NLSP decay, and 
y NLSP = 77 N LSp/n5 G is the NLSP number density before they decay, normalized to 
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the number density of background photons n® G = 2((3)T 3 /n 2 . We define hadronic 
energy release similarly as Chad = £had*NLSP- 

NLSP velocities are negligible when they decay, and so the potentially visible en- 
ergy is 

E s ee mgjLSP ~ m ^ . (53) 

2mNLSP 

For the photino case, S = 7. At leading order, all of the initial photon energy is 
deposited in an electromagnetic shower, and so 

E E M = £7 , £had ^ . (54) 

For the stau case, 

^EM ~ 7^E T — E T , £ had = , (55) 

where the range in £em results from the possible variation in electromagnetic energy 
from ix^~ and v decay products. The precise value of £em is in principle calculable once 
the stau's chirality and mass, and the superWIMP mass, are specified. However, as the 
possible variation in e EM is not great relative to other effects, we will simply present 
results below for the representative value of e EM = \E T . 

The lifetimes and energy releases in the photino and stau NLSP scenarios are given 
in Fig. f° r a range of (m NL sp, Am). For natural weak-scale values of these param- 
eters, the lifetimes and energy releases in the neutralino and stau scenarios are similar, 
with lifetimes of about a year, in accord with the rough estimate of Eq. (l46b . and energy 
releases of 

Cem ~ 10- 9 GeV . (56) 
Such values have testable implications, as we now discuss. 



4.5.2 Big Bang Nucleosynthesis 

Big Bang nucleosynthesis predicts primordial light element abundances in terms of one 
free parameter, the baryon-to-photon ratio r\ = n B /n 7 . At present, the observed D, 
4 He, 3 He, and 7 Li abundances may be accommodated for baryon-to-photon ratios in 
the range 132 

7710 = 77/Kr 10 = 2.6 - 6.2 . (57) 

(See Fig. |23l) In light of the difficulty of making precise theoretical predictions and 
reducing (or even estimating) systematic uncertainties in the observations, this consis- 
tency is a well-known triumph of standard Big Bang cosmology. 
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10 6 10 9 10 12 10 15 10 18 10 6 10 9 10 12 10 15 10 18 
t wimp (sec) t wimp (sec) 



Fig. 22. Predicted values of NLSP lifetime r and electromagnetic energy release £em = 
^em^nlsp in the 7 (left) and f (right) NLSP scenarios for = 1 GeV, 10 GeV, 
100 TeV (top to bottom) and Am = m NLSP - m & = 1 TeV, 100 GeV, . . . , 100 MeV 
(left to right). For the f NLSP scenario, we assume £em — \E T . From Ref. ll09[ 

At the same time, given recent and expected advances in precision cosmology, the 
standard BBN picture merits close scrutiny. Recently, BBN baryometry has been sup- 
plemented by CMB data, which alone yields i]i — 6.1 ± 0.4. 1 Observations of deu- 
terium absorption features in spectra from high redshift quasars imply a primordial D 
fraction of D/H = 2.781^ x lO" 5 .^ Combined with standard BBN calculations,^ 
this yields 7]i = 5.9 ± 0.5. The remarkable agreement between CMB and D baryome- 
ters has two new implications for scenarios with late-decaying particles. First, assuming 
there is no fine-tuned cancellation of unrelated effects, it prohibits significant entropy 
production between the times of BBN and decoupling. Second, the CMB measurement 
supports determinations of i] from D, already considered by many to be the most re- 
liable BBN baryometer. It suggests that if D and another BBN baryometer disagree, 
the "problem" lies with the other light element abundance — either its systematic un- 
certainties have been underestimated, or its value is modified by new astrophysics or 
particle physics. At present BBN predicts a 7 Li abundance significantly greater ob- 
served. This disagreement may therefore provide specific evidence for late-decaying 
particles in general, and gravitino dark matter in particular. 

Given the overall success of BBN, the first implication for new physics is that it 
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Fig. 23. Bounds on the baryon density f2#/i 2 from BBN (left, from Ref. 11321) and the 
CMB (right, from Ref. 11331) . The new and extremely precise CMB constraint favors 
the BBN VtBh 2 determination from deuterium and implies that the 7 Li abundance is 
anomalously low. 

should not drastically alter any of the light element abundances. This requirement re- 
stricts the amount of energy released at various times in the history of the universe. A 
recent analysis of electromagnetic cascades finds that the shaded regions of Fig.l24lare 
excluded by such considerations. 127 The various regions are disfavored by the follow- 
ing conservative criteria: 

Dlow : D/H < 1.3 x 1(T 5 (58) 
D high : D/H > 5.3 x 1(T 5 (59) 
4 He low : Y p < 0.227 (60) 
7 Li low : 7 Li/H < 0.9 x 10~ 10 . (61) 



A subset of superWIMP predictions from Fig.[22]is superimposed on this plot. The 
subset is for weak-scale and Am, the most natural values, given the independent 
motivations for new physics at the weak scale. The BBN constraint eliminates some of 
the region predicted by the superWIMP scenario, but regions with m-NLSP, m c ~ -^weak 



45 




Fig. 24. The grid gives predicted values of NLSP lifetime r and electromagnetic 
energy release Cem = ^em^nlsp in the 7 (left) and f (right) NLSP scenarios for 
m<5 = 100 GeV, 300 GeV, 500 GeV, 1 TeV, and 3 TeV (top to bottom) and 
Am = m NLSP -m d = 600 GeV, 400 GeV, 200 GeV, and 100 GeV (left to right). 
For the f NLSP scenario, we assume e EM = \E T . BBN constraints exclude the shaded 
regions. 123 The best fit region with (r, Cem) ~ (3 x 10 6 s, 10" 9 GeV), where 7 Li is 
reduced to observed levels by late decays of NLSPs to gravitinos, is given by the circle. 
FromRef.fTD^l 

remain viable. 

The 7 Li anomaly discussed above may be taken as evidence for new physics, how- 
ever. To improve the agreement of observations and BBN predictions, it is necessary to 
destroy 7 Li without harming the concordance between CMB and other BBN determi- 
nations of f]. This may be accomplished for (r, (em) ~ (3 x 10 6 s, 10~ 9 GeV).^ This 
"best fit" point is marked in Fig. EH The amount of energy release is determined by the 
requirement that 7 Li be reduced to observed levels without being completely destroyed 
- one cannot therefore be too far from the " 7 Li low" region. In addition, one cannot 
destroy or create too much of the other elements. 4 He, with a binding threshold energy 
of 19.8 MeV, much higher than Lithium's 2.5 MeV, is not significantly destroyed. On 
the other hand, D is loosely bound, with a binding energy of 2.2 MeV. The two primary 
reactions are D destruction through 7D — > np and D creation through 7 4 He — > DD. 
These are balanced in the channel of Fig.|24|between the "low D" and "high D" regions, 
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and the requirement that the electromagnetic energy that destroys Li not disturb the D 
abundance specifies the preferred decay time r ~ 3 x 10 6 s. 

Without theoretical guidance, this scenario for resolving the 7 Li abundance is rather 
fine-tuned: possible decay times and energy releases span tens of orders of magnitude, 
and there is no motivation for the specific range of parameters required to resolve BBN 
discrepancies. In the superWIMP scenario, however, both r and (em are specified: the 
decay time is necessarily that of a gravitational decay of a weak-scale mass particle, 
leading to Eq. d4oT) . and the energy release is determined by the requirement that su- 
perWIMPs be the dark matter, leading to Eq. <l56b . Remarkably, these values coincide 
with the best fit values for r and (em- More quantitatively, we note that the grids of 
predictions for the 7 and f scenarios given in Fig.|24|cover the best fit region. Current 
discrepancies in BBN light element abundances may therefore be naturally explained 
by gravitino dark matter. 

This tentative evidence may be reinforced or disfavored in a number of ways. 
Improvements in the BBN observations discussed above may show if the 7 Li abun- 
dance is truly below predictions. In addition, measurements of 6 Li/H and 6 Li/ 7 Li 
may constrain astrophysical depletion of 7 Li and may also provide additional evi- 
dence for late decaying particles in the best fit re gi on- EHEMI13n3S3 Finally, if the 
best fit region is indeed realized by NLSP — > G decays, there are a number of other 
testable implications for cosmology and particle physics. We discuss one of these in 
the following section. Additional discussion, including diffuse photon signals, the im- 
plications of hadronic energy release, and novel collider analyses, may be found in 
Ref s . 110811091138113911401141114211431 

4.5.3 The Cosmic Microwave Background 

The injection of electromagnetic energy may also distort the frequency dependence of 
the CMB black body radiation. For the decay times of interest, with redshifts z ~ 10 5 — 
10 7 , the resulting photons interact efficiently through 7e~ — > 7e~, but photon number is 
conserved, since double Compton scattering — > 77c - and thermal bremsstrahlung 
eX — > eX'j, where X is an ion, are inefficient. The spectrum therefore relaxes to 
statistical but not thermodynamic equilibrium, resulting in a Bose-Einstein distribution 
function 

= eE/{k T )+ , _ l . (62) 

with chemical potential /U ^ 0. 
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Fig. 25. Contours of fj,, parameterizing the distortion of the CMB from a Planckian 
spectrum, in the (r, £em) plane. Regions predicted by the gravitino dark matter sce- 
nario, and BBN excluded and best fit regions are given as in Fig. |24j From Ref. 11091 



For the low values of baryon density currently favored, the effects of double Comp- 
ton scattering are more significant than those of thermal bremsstrahlung. The value of 
the chemical potential [i may therefore be approximated for small energy releases by 
the analytic expressiorP^ 



where 



H = 8.0 x 10" 



r dC = 6.1 x 10 6 s 



T 



10 6 s 



T 
2.725 K 



EAl 



10- 9 GeV 



Q R h 2 



0.022 



-(W^) 5/4 



2 V 



0.88 



(63) 



(64) 



In Fig.|25]we show contours of chemical potential The current bound is /i < 
9 x 10 -5 We see that, although there are at present no indications of deviations 

from black body, current limits are already sensitive to the superWIMP scenario, and 
particularly to regions favored by the BBN considerations described in Sec. 14.5.21 In 
the future, the Diffuse Microwave Emission Survey (DIMES) may improve sensitivities 
to /j w 2 x io _6 P2l DIMES will therefore probe further into superWIMP parameter 
space, and will effectively probe all of the favored region where the 7 Li underabundance 
is explained by gravitino dark matter. 
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4.6 Summary 

• The gravitino mass is determined by the scale of supersymmetry breaking and 
may be anywhere in the range from eV to TeV. In supergravity theories, its mass 
is at the weak scale and its couplings are suppressed by the Planck scale, and so 
extremely weak. 

• If gravitinos are produced as a thermal relic, their mass is bounded by overclosure 
to be vtlq < keV if they are stable, and by BBN to be > 10 TeV if they are 
unstable. 

• Gravitinos may be produced after inflation during reheating. For stable weak- 
scale gravitinos, overclosure places an upper bound on the reheat temperature of 
the order of 10 10 GeV. 

• Weak-scale gravitinos may also be produced in NLSP decays at time t ~ 10 4 — 
10 8 s. In this case, gravitinos may be dark matter, naturally inheriting the desired 
relic density. Gravitino dark matter is undetectable by conventional direct and 
indirect dark matter searches, but may be discovered through its imprint on early 
universe signals, such as BBN and the CMB. 

5 Prospects 

We have now discussed a wide variety of cosmological implications of supersymmetry. 
If discoveries are made in astrophysical and cosmological observations, what are the 
prospects for determining if this new physics is supersymmetry? Put more generally, 
what are the prospects for a microscopic understanding of the dark universe? Such 
questions are grand, and their answers speculative. Nevertheless, some lessons may 
be drawn even now. As we will see, even in the best of cases, we will need diverse 
experiments from both particle physics and cosmology to explore this frontier. 

5.1 The Particle Physics/Cosmology Interface 

As a case study, let us confine our discussion to one topic: neutralino dark matter. We 
assume that non-baryonic dark matter is in fact neutralinos. If this is so, what are the 
prospects for establishing this, and what tools will we need? 

It is first important to recognize the limitations of both cosmology and particle 
physics when taken separately: 



49 



• Cosmological observations and astrophysical experiments cannot discover super- 
symmetry. As noted in Sec. \\\ cosmological data leaves the properties of dark 
matter largely unconstrained. If dark matter is discovered in direct or indirect de- 
tection experiments, its mass and interaction strengths will be bounded but only 
very roughly at first. (For example, the region favored by the DAMA signal spans 
factors of a few in both mass and interaction strength; see Fig. [14|) These con- 
straints will be sharpened by follow-up experiments. However, the microscopic 
implications of such experiments are clouded by significant astrophysical ambi- 
guities, such as the dark matter velocity distribution, halo profiles, etc. Even with 
signals in a variety of direct and indirect detection experiments, it is unlikely that 
dark matter properties will be constrained enough to differentiate supersymmetry 
from other reasonable possibilities. 

• Particle experiments cannot discover dark matter. If weak- scale superpartners 
exist, particle colliders will almost certainly be able to discover at least some 
of them. However, even if they find all of them, the dark matter candidate will 
most likely appear only as missing energy and momentum. Furthermore, collider 
experiments can only test the stability of such particles up to lifetimes of ~ 1CT 7 s. 
As we have seen in Sec.HJ lifetimes of a year or more are perfectly natural in well- 
motivated models of new physics. The conclusion that a particle seen in collider 
experiments is the dark matter therefore requires an unjustified extrapolation of 
24 orders of magnitude in the particle's lifetime. 

Through the combination of both approaches, however, it is possible that a cohe- 
sive and compelling theory of dark matter will emerge. A schematic picture of the 
combined investigation of neutralino dark matter is given in Fig.|26l 147 Working from 
the bottom, cosmological observations have already determined the relic density with 
some precision. Future observations, such as by the Planck satellite]^ are likely to 
reduce uncertainties in the relic density determination to the 1 % level, given now stan- 
dard cosmological assumptions. Astrophysical experiments may also detect dark matter 
either directly through its interactions with ordinary matter or indirectly through its an- 
nihilation decay products. Such data, combined with astrophysical inputs such as the 
dark matter halo profile and local density, will provide information about the strength 
of xiV scattering and XX annihilation. 

At the same time, working from the top of Fig. |26l colliders will discover super- 
symmetry and begin to determine the parameters of the weak-scale Lagrangian. These 
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Fig. 26. The road to a microscopic understanding of neutralino dark matter. 

parameters will, in principle, fix the neutralino 's thermal relic density, the %iV scatter- 
ing cross section, and the neutralino pair annihilation rates. Completion of this program 
at a high level of precision, followed by detailed comparison with the measured relic 
density and detection rates from cosmology and astrophysics will provide a great deal 
of information about the suitability of neutralinos as dark matter candidates. 



Clearly data from particle colliders will be required to identify neutralino dark matter. 
The requirements for colliders depend sensitively on what scenario is realized in nature. 
As examples, consider some of the cosmologically preferred regions of minimal super- 
gravity discussed in Sec. 13.21 In the bulk region, one must verify that the neutralino 
is Bino-like and must determine the masses of sfermions that appear in the t-channel 
annihilation diagrams. In the focus point region, the neutralino's gaugino-ness must be 
precisely measured, whereas in the A funnel region, a high precision measurement of 
rriA — 2m x is critical. Finally, for the co-annihilation region, there is extreme sensitivity 
to the f~x mass splitting. Measurements below the GeV level are required to accurately 
determine the predicted thermal relic density. 

Let us consider the bulk region scenario in more detail. Not all sfermion masses 
need be measured. For example, if the right-handed sleptons are light, they typically 



5.2 The Role of Colliders 
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Fig. 27. The ratio of the true iluMh 2 to that calculated with the Ir i-channel diagrams 
in the (m, M) plane, where m and M are the universal scalar and gaugino masses of 
minimal supergravity, respectively. From Ref . 11491 

give the dominant contribution, since these have the largest hypercharge Y and the anni- 
hilation diagram is proportional to F 4 . In such cases, measurements of m~ [R and lower 
bounds on left-handed slepton and squark masses will provide a reasonable starting 
point. 

The possibility of doing this at the LHC has been considered in Ref. 11491 In much 
of the bulk region, the cascade decay q L — > xll ~ * W<7 ~~ * l + l~XxQ 1S open. Kinematic 
endpoints may then be used to determine the Ir and x\ masses precisely. Assuming that 
the lightest neutralino is Bino-like, one may then estimate the relic density, keeping only 
Ir exchange diagrams. As shown in Fig.|27J this provides an estimate accurate to about 
~ 20% in minimal supergravity. Following this, one would then need to determine 
the gaugino-ness of the lightest neutralino and set lower bounds on the other sfermion 
masses. 

At a linear collider, one may establish that the new particles being produced are 
supersymmetric by measuring their dimensionless couplings. One may then go on to 
determine the gaugino-ness of the LSP in a model-independent manner. For exam- 
ple, the cross section a{e + e r R — > x + X ) nearly vanishes for gaugino-like charginos. It 
therefore provides a sensitive measure of chargino gaugino-ness. (See Fig.^El) Com- 
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Fig. 28. Contours of constant cross section a(e + e R — > x + X ) f° r a V$ — 500 GeV 
linear collider. From Ref . 11501 

bined with kinematic measurements of the chargino mass, the parameters M 2 and /i 
may be measured precisely. Further measurements can use these results to pinpoint Mi 
and tan (3, and thereby the gaugino-ness of the LSP. Precisions of ~ 1% or better are 
possible, translating into predictions for relic densities and dark matter cross sections 
that will match the precision expected from cosmological data. 

5.3 Synthesis 

If the relic density and interaction strengths as determined by astrophysics and cosmol- 
ogy agree with the predictions of particle physics with high precision, this agreement 
will provide strong evidence that the dark matter is in fact supersymmetric. It will im- 
ply that we understand the history of the universe back to the freeze out temperature 
~ 10 GeV, or times t ~ 10~ 8 s. Recall that our current knowledge of the history of 
the universe is on sure footing only back to Big Bang nucleosynthesis at temperatures 
of ~ MeV, or times t ~ 1 s. Dark matter studies could therefore provide the necessary 
evidence to push back our knowledge of the universe another 8 orders of magnitude in 
time, a formidable achievement. 

On the other hand, the determinations of relic density and dark matter interaction 
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strengths by particle physics and cosmology may not agree. Progress then has many 
possible paths. If the disparity is great, one might look to other dark matter candidates, 
such as the axiori 151 1 15 -^ 1 ^ 1 or other supersymmetric possibilities! 154 ' 155 1 156 1 110 1 111 1 11 ^ If 
the relic density determinations are reasonably close, one might explore the possibility 
that the neutralino is not stable, but deposits much of its relic density in a gravitino LSP, 
as discussed in Sec. 14.41 

Alternatively, one might look to non-standard cosmologies for a resolution. The 
identification of the thermal relic density with the present day cold dark matter density 
is subject to cosmological assumptions. The calculation of the thermal relic density as- 
sumes that the dominant source of dark matter is from dark matter particles falling out 
of thermal equilibrium. It is possible, however, that the bulk of the dark matter is cre- 
ated not through thermal equilibrium and freeze out, but through the out-of-equilibrium 
decay of a supermassive particle. The actual relic density would then be greater than the 
thermal relic density. The thermal relic density calculation also assumes that nothing 
unusual happens once the dark matter is produced at temperatures of T ~ (9(10) GeV. 
Large entropy production by late-decaying particles may dilute calculated relic den- 
sities. In this case, the actual relic density would be less than the naive thermal relic 
density. The bottom line is that the cold dark matter density obtained following the path 
from the bottom of Fig.|26]need not coincide with the thermal relic density obtained by 
following the path from the top. Instead, discrepancies might provide new insights into 
the history of our universe. 

In a similar vein, the neutralino-nucleon cross sections need not match the dark 
matter detection rates. As stressed above, this correspondence requires astrophysical 
assumptions. The uncertainties and problems associated with these issues have been 
discussed extensively! 28 ' 29 1 31 ^ It is possible, however, that the relic densities, as deter- 
mined independently by particle physics and cosmology, agree to 1%, but the detection 
rates differ. One would then be confident that neutralinos are the dark matter and parti- 
cle physics uncertainties would be eliminated, allowing detection experiments to probe 
astrophysics. For example, direct detection rates would then provide information about 
the local dark matter density and velocity distributions, and indirect detection rates 
would provide information about halo profiles. The synergy between cosmology and 
particle physics would then truly come full circle. 
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5.4 Summary 

A microscopic understanding of the dark universe is a challenging goal. As an exam- 
ple, we have focused here on prospects for a fundamental description of dark matter. 
Cosmological measurements, although able to bound total energy densities, cannot tell 
us much about the dark matter's microscopic properties. On the other hand, particle 
physics experiments may produce dark matter and may measure its properties rather 
precisely, but cannot never establish its stability on cosmological time scales. It is only 
through the combination of approaches in particle physics, astrophysics, and cosmol- 
ogy that the identity of dark matter will be uncovered. The task requires many diverse 
experiments, and will likely take decades to complete. Nevertheless, if any of the con- 
nections between the weak scale and cosmology described here are realized in nature, 
one would be hard-pressed to envision a more exciting era of discovery than the coming 
years. 

6 Acknowledgments 

I am grateful to the organizers of the 2003 SLAC Summer Institute, where this ma- 
terial was first presented, and to the students and participants for their enthusiasm and 
lively discussions. I thank the many colleagues, especially Konstantin Matchev, Arvind 
Rajaraman, Fumihiro Takayama, and Frank Wilczek, who have helped to shape my un- 
derstanding of the topics presented here. This work was supported in part by National 
Science Foundation CAREER Award PHY-0239817 and in part by the Alfred P. Sloan 
Foundation. 

References 

[1] D. N. Spergel et ah, "First Year Wilkinson Microwave Anisotropy 
Probe (WMAP) Observations: Determination of Cosmological Parameters," 
astrcTph/0302209| 

[2] M. Tegmark et al. [SDSS Collaboration], "Cosmological parameters from SDSS 
and WMAP," |astro-ph/03 T0723l 

[3] For proper introductions to supersymmetry, see, for example, H. J. Muller- 
Kirsten and A. Wiedemann, Supersymmetry: An Introduction With Concep- 
tual And Calculation^ Details, Print-86-0955 (KAIS ERS LAUTERN) ; J. Wess 



55 



and J. Bagger, Super symmetry and super gravity, 2nd edition (Princeton Uni- 
versity Press, Princeton, NT, 1992); H. E. Haber, "Introductory low-energy 
supersymmetry," hep-ph/9306207 ; X. Tata, "Super symmetry: Where it is 
and how to find it," hep-ph/95 10287, M. Drees, "An introduction to super- 
symmetry," hep-ph/961 1409; J. D. Lykken, "Introduction to supersymmetry," 
hep-th/96 12114} S. P. Martin, "A supersymmetry primer," hep-ph/9709356; 
S. Weinberg, The Quantum Theory Of Fields. Vol. 3: Supersymmetry (Cam- 
bridge University Press, Cambridge, UK, 2000); N. Polonsky, Supersymmetry: 
Structure and phenomena. Extensions of the standard model, Lect. Notes Phys. 
M68, 1 (2001) | |hep^p h/0108236|. 

[4] R. Haag, J. T. Lopuszanski and M. Sohnius, "All Possible Generators Of Super- 
symmetries Of The S Matrix," Nucl. Phys. B 88, 257 (1975). 

[5] S. Dimopoulos, S. Raby and F. Wilczek, "Supersymmetry And The Scale Of 
Unification," Phys. Rev. D 24, 1681 (1981). 

[6] A. H. Chamseddine, R. Arnowitt and P. Nath, "Locally Supersymmetric Grand 
Unification," Phys. Rev. Lett. 49, 970 (1982). 

[7] R. Barbieri, S. Ferrara and C. A. Savoy, "Gauge Models With Spontaneously 
Broken Local Supersymmetry," Phys. Lett. B 119, 343 (1982). 

[8] N. Ohta, "Grand Unified Theories Based On Local Supersymmetry," Prog. 
Theor. Phys. 70, 542 (1983). 

[9] L. J. Hall, J. Lykken and S. Weinberg, "Supergravity As The Messenger Of Su- 
persymmetry Breaking," Phys. Rev. D 27 (1983) 2359. 

[10] L. Alvarez-Gaume, J. Polchinski and M. B. Wise, "Minimal Low-Energy Super- 
gravity," Nucl. Phys. B 221, 495 (1983). 

[11] K. A. Olive, "TASI lectures on dark matter," |astro~- ph/030 1 505 

[12] J. L. Feng, K. T. Matchev and F. Wilczek, "Neutralino dark matter in focus point 
supersymmetry," Phys. Lett. B 482, 388 (2000) [hep-ph/0004043|. 

[13] H. Goldberg, "Constraint On The Photino Mass From Cosmology," Phys. Rev. 
Lett. 50, 1419 (1983). 

[14] J. R. Ellis, J. S. Hagelin, D. V. Nanopoulos, K. A. Olive and M. Srednicki, "Su- 
persymmetric Relics From The Big Bang," Nucl. Phys. B 238, 453 (1984). 



56 



[15] For proper introductions to dark matter, see, for example, Refs. 1 1 61 1711 1 and 
L. Bergstrom, "Non-baryonic dark matter: Observational evidence and detection 
methods," Rept. Prog. Phys. 63, 793 (2000) |hep-ph/0002126|. 

[16] G. Jungman, M. Kamionkowski and K. Griest, "Super symmetric dark matter," 
Phys. Rept. 267, 195 (1996) [hep-ph/9506380|. 

[17] E. W. Kolb and M. S. Turner, The Early Universe (Addison-Wesley, Redwood 
City, CA, 1990). 

[18] M. Drees and M. M. Nojiri, "The Neutralino relic density in minimal N=l su- 
pergravity," Phys. Rev. D 47, 376 (1993) [hep-ph/9207234]. 

[19] P. Gondolo, J. Edsjo, L. Bergstrom, P. Ullio and E. A. Baltz, "DarkSUSY: A nu- 
merical package for dark matter calculations in the MSSM," astro-ph/00 12234 

[20] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, "Supersymmetric dark 
matter in light of WMAP," Phys. Lett. B 565, 176 (2003) |hep-ph/0303043|. 

[21] J. L. Feng, K. T. Matchev and T. Moroi, "Multi-TeV scalars are natural in mini- 
mal supergravity," Phys. Rev. Lett. 84, 2322 (2000) |hep-ph/9908309|. 

[22] J. L. Feng, K. T. Matchev and T. Moroi, "Focus points and naturalness in super- 
symmetry," Phys. Rev. D 61, 075005 (2000) |hep-ph/9909334|. 

[23] J. L. Feng and K. T. Matchev, "Focus point supersymmetry: Proton decay, flavor 
and CP violation, and the Higgs boson mass," Phys. Rev. D 63, 095003 (2001) 
[hep-ph/ 0011356| l. 

[24] H. Baer, C. Balazs and A. Belyaev, "Neutralino relic density in minimal super- 
gravity with co-annihilations," JHEP 0203, 042 (2002) [hep-ph/0202076|. 

[25] H. Baer, C. Balazs and A. Belyaev, "Relic density of neutralinos in minimal 
supergravity," hep-ph/02 11213 

[26] P. Binetruy, G. Girardi and P. Salati, "Constraints On A System Of Two Neutral 
Fermions From Cosmology," Nucl. Phys. B 237, 285 (1984). 

[27] K. Griest and D. Seckel, "Three Exceptions In The Calculation Of Relic Abun- 
dances," Phys. Rev. D 43, 3191 (1991). 

[28] P. Sikivie, I. I. Tkachev and Y. Wang, "The secondary infall model of galactic 
halo formation and the spectrum of cold dark matter particles on earth," Phys. 
Rev. D 56, 1863 (1997) |astro-ph/9609022|. 



57 



[29] M. Brhlik and L. Roszkowski, "WIMP velocity impact on direct dark matter 
searches," Phys. Lett. B 464, 303 (1999) [hep-ph/9903468|. 

[30] P. Belli, R. Bernabei, A. Bottino, F. Donato, N. Fornengo, D. Prosperi and 
S. Scopel, "Extending the DAMA annual-modulation region by inclusion of 
the uncertainties in astrophysical velocities," Phys. Rev. D 61, 023512 (2000) 
[hep-ph/99033UTJ. 

[31] G. Gelmini and P. Gondolo, "WIMP annual modulation with opposite phase in 
late-infall halo models," Phys. Rev. D 64, 023504 (2001) [hep-ph/ 00T23T5l . 

[32] C. Calcaneo-Roldan and B. Moore, "The surface brightness of dark matter: 
Unique signatures of neutralino annihilation in the galactic halo," Phys. Rev. 
D 62, 123005 (2000) [astro-ph/00 T0036| |. 

[33] M. W. Goodman and E. Witten, "Detectability Of Certain Dark-Matter Candi- 
dates," Phys. Rev. D 31, 3059 (1985). 

[34] J. R. Ellis, A. Ferstl and K. A. Olive, "Re-evaluation of the elastic scattering of 
supersymmetric dark matter," Phys. Lett. B 481, 304 (2000) [hep-ph/0001005]. 

[35] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, "Remarks On Higgs - Boson 
Interactions With Nucleons," Phys. Lett. B 78, 443 (1978). 

[36] H. Baer, C. Balazs, A. Belyaev and J. O'Farrill, "Direct detection of dark matter 
in supersymmetric models," JCAP 0309, 007 (2003) |hep-ph/0305r9"T1 |. 

[37] R. Bernabei et al. [DAMA Collaboration], "Search for WIMP annual modula- 
tion signature: Results from DAMA / NaI-3 and DAMA / NaI-4 and the global 
combined analysis," Phys. Lett. B 480, 23 (2000). 

[38] A. Benoit et al, "Improved exclusion limits from the EDELWEISS WIMP 
search," Phys. Lett. B 545, 43 (2002) [astro-ph/0206271]. 

[39] D. S. Akerib et al. [CDMS Collaboration], "New results from the Cryo- 
genic Dark Matter Search experiment," Phys. Rev. D 68, 082002 (2003) 
|hep-ex/030600f). 

[40] D. S. Akerib et al. [CDMS Collaboration], "First results from the cryogenic dark 
matter search in the Soudan Underground Lab," astro-ph/0405033. 

[41] C. J. Copi and L. M. Krauss, "Comparing WIMP interaction rate detec- 
tors with annual modulation detectors," Phys. Rev. D 67, 103507 (2003) 
[astro-ph/0208010]. 



58 



[42] K. Freese, P. Gondolo and H. J. Newberg, "Detectability of weakly interacting 
massive particles in the Sagittarius dwarf tidal stream," astro-ph/0309279. 

[43] P. Ullio, M. Kamionkowski and P. Vogel, "Spin dependent WEVIPs in DAMA?," 
JHEP 0107, 044 (2001) [hep-ph/0010036|. 

[44] D. R. Smith and N. Weiner, "Inelastic dark matter," Phys. Rev. D 64, 043502 
(2001) ]hep-ph/0101O31 |. 

[45] A. Kurylov and M. Kamionkowski, "Generalized analysis of weakly-interacting 
massive particle searches," Phys. Rev. D 69, 063503 (2004) |Kej>jph/0307T55|. 

[46] D. Tucker-Smith and N. Weiner, "The status of inelastic dark matter," 
|hep-ph/0402065| 

[47] S. Rudaz and F. W. Stecker, "Cosmic Ray Anti-Protons, Positrons And Gamma- 
Rays From Halo Dark Matter Annihilation," Astrophys. J. 325, 16 (1988). 

[48] A. J. Tylka, "Cosmic Ray Positrons From Annihilation Of Weakly Interacting 
Massive Particles In The Galaxy," Phys. Rev. Lett. 63, 840 (1989). 

[49] M. S. Turner and F. Wilczek, "Positron Line Radiation From Halo Wimp Anni- 
hilations As A Dark Matter Signature," Phys. Rev. D42, 1001 (1990). 

[50] M. Kamionkowski and M. S. Turner, "A Distinctive positron feature from heavy 
WIMP annihilations in the galactic halo," Phys. Rev. D43, 1774 (1991). 

[51] E. A. Baltz and J. Edsjo, "Positron propagation and fluxes from neutralino anni- 
hilation in the halo," Phys. Rev. D59, 02351 1 (1999) [astro-ph/9808243|. 

[52] I. V. Moskalenko and A. W. Strong, "Positrons from particle dark-matter anni- 
hilation in the galactic halo: Propagation Green's functions," Phys. Rev. D60, 
063003 (1999) [astro-ph/9905283]. 

[53] F. W. Stecker and A. J. Tylka, "The Cosmic Ray Anti-Proton Spectrum From 
Dark Matter Annihilation And Its Astrophysical Implications: A New Look," 
Ap. J. 336, L51 (1989). 

[54] P. Chardonnet, G. Mignola, P. Salati and R. Taillet, "Galactic diffusion and the 
antiproton signal of supersymmetric dark matter," Phys. Lett. B384, 161 (1996) 
[astro-ph/9606174]. 

[55] L. Bergstrom, J. Edsjo and P. Ullio, "Cosmic antiprotons as a probe for super- 
symmetric dark matter?," astro-ph/9902012. 



59 



[56] J. W. Bieber, R. A. Burger, R. Engel, T. K. Gaisser, S. Roesler and 
T. Stanev, "Antiprotons at solar maximum," Phys. Rev. Lett. 83, 674 (1999) 
[Jastro-ph/9903 163 1 . 

[57] F. Donato, N. Fornengo and R Salati, "Antideuterons as a signature of supersym- 
metric dark matter," Phys. Rev. D 62, 043003 (2000) |hep-ph/9904481 1. 

[58] G. Servant and T. M. P. Tait, "Is the lightest Kaluza-Klein particle a viable dark 
matter candidate?" Nucl. Phys. B 650, 391 (2003) |hep-ph/0206071|. 

[59] H. C. Cheng, J. L. Feng and K. T. Matchev, "Kaluza-Klein dark matter," Phys. 
Rev. Lett. 89, 211301 (2002) [ |hep-ph/0207T23J . 

[60] R. J. Protheroe, "On The Nature Of The Cosmic Ray Positron Spectrum," Astro- 
phys. J. 254, 391 (1982). 

[61] S. W. Barwick et al. [HEAT Collaboration], "Cosmic ray positrons at 
high-energies: A New measurement," Phys. Rev. Lett. 75, 390 (1995) 
|astro-ph/9505141| |. 

[62] S. W. Barwick et al. [HEAT Collaboration], "Measurements of the cosmic- 
ray positron fraction from 1-GeV to 50-GeV," Astrophys. J. 482, L191 (1997) 
|astro-ph/9703192| . 

[63] S. Coutu et al. [HEAT-pbar Collaboration], "Positron Measurements with the 
HEAT-pbar Instrument," in Proceedings of the 27th International Cosmic Ray 
Conference (2001). 

[64] G. L. Kane, L. T. Wang and J. D. Wells, "Super symmetry and the positron excess 
in cosmic rays," Phys. Rev. D 65, 057701 (2002) [hep-ph/ 0l08T38| . 

[65] E. A. Baltz, J. Edsjo, K. Freese and P. Gondolo, "The cosmic ray positron excess 
and neutralino dark matter," Phys. Rev. D 65, 0635 1 1 (2002) [ astro-ph/01093 1 8 1 . 

[66] G. L. Kane, L. T. Wang and T. T. Wang, "Supersymmetry and the cosmic ray 
positron excess," Phys. Lett. B 536, 263 (2002) [hep-ph/0202156]. 

[67] D. Hooper, J. E. Taylor and J. Silk, "Can supersymmetry naturally explain the 
positron excess?," hep-ph/03 12076, 

[68] F. W. Stecker, "Gamma-Ray Constraints On Dark Matter Reconsidered," Phys. 
Lett. B 201, 529 (1988). 

[69] F. W. Stecker and A. J. Tylka, "Spectra, Fluxes And Observability Of Gamma- 
Rays From Dark Matter Annihilation In The Galaxy," Astrophys. J. 343, 169 
(1989). 



60 



[70] L. Bergstrom, P. Ullio and J. H. Buckley, "Observability of gamma rays from 
dark matter neutralino annihilations in the Milky Way halo," Astropart. Phys. 9, 
137 (1998) |astro-ph/97l25Tq . 

[71] V. S. Berezinsky, A. V. Gurevich and K. P. Zybin, "Distribution of dark matter 
in the galaxy and the lower limits for the masses of supersymmetric particles," 
Phys. Lett. B294, 221 (1992). 

[72] V. Berezinsky, A. Bottino and G. Mignola, "High-energy gamma radiation from 
the galactic center due to neutralino annihilation," Phys. Lett. B325, 136 (1994) 
[hep-ph/ 9402215| . 

[73] M. Urban, A. Bouquet, B. Degrange, P. Fleury, J. Kaplan, A. L. Melchior and 
E. Pare, "Searching for TeV dark matter by atmospheric Cerenkov techniques," 
Phys. Lett. B293, 149 (1992) [hep-ph/9208255|. 

[74] P. Gondolo and J. Silk, "Dark matter annihilation at the galactic center," Phys. 
Rev. Lett. 83, 1719 (1999) | |astro-ph/ 990639il- 

[75] L. Bergstrom, J. Edsjo and P. Ullio, "Possible Indications of a Clumpy Dark 
Matter Halo," Phys. Rev. D58, 083507 (1998) [astro-ph/9804050|. 

[76] L. Bergstrom, J. Edsjo, P. Gondolo and P. Ullio, "Clumpy neutralino dark mat- 
ter," Phys. Rev. D59, 043506 (1999) (astn>ph/9806072|. 

[77] E. A. Baltz, C. Briot, P. Salati, R. Taillet and J. Silk, "Detection of neutralino 
annihilation photons from external galaxies," Phys. Rev. D61, 023514 (2000) 
|astro-ph/9909112| |. 

[78] S. Rudaz and F. W. Stecker, "On The Observability Of The Gamma-Ray Line 
Flux From Dark Matter Annihilation," Astrophys. J. 368, 406 (1991). 

[79] L. Bergstrom and P. Ullio, "Full one-loop calculation of neutralino annihilation 
into two photons," Nucl. Phys. B504, 27 (1997) [hep-ph/97U(>232). 

[80] Z. Bern, P. Gondolo and M. Perelstein, "Neutralino annihilation into two pho- 
tons," Phys. Lett. B411, 86 (1997) [hep-ph/9706538|. 

[81] P. Ullio and L. Bergstrom, "Neutralino annihilation into a photon and a Z boson," 
Phys. Rev. D57, 1962 (1998) |hep-ph/9707333|. 

[82] V. S. Berezinsky, A. Bottino and V. de Alfaro, "Is it possible to detect the gamma 
ray line from neutralino-neutralino annihilation?" Phys. Lett. B274, 122 (1992). 



61 



[83] K. Freese, "Can Scalar Neutrinos Or Massive Dirac Neutrinos Be The Missing 
Mass?," Phys. Lett. B167, 295 (1986). 

[84] L. M. Krauss, M. Srednicki and F. Wilczek, "Solar System Constraints And Sig- 
natures For Dark Matter Candidates," Phys. Rev. D33, 2079 (1986). 

[85] T. K. Gaisser, G. Steigman and S. Tilav, "Limits On Cold Dark Matter Candi- 
dates From Deep Underground Detectors," Phys. Rev. D34, 2206 (1986). 

[86] A. Gould, J. A. Frieman and K. Freese, "Probing The Earth With Wimps," Phys. 
Rev. D39, 1029 (1989). 

[87] A. Bottino, N. Fornengo, G. Mignola and L. Moscoso, "Signals of neutralino 
dark matter from earth and sun," Astropart. Phys. 3, 65 (1995) [hep-ph/9408391 1. 

[88] V. Berezinsky, A. Bottino, J. R. Ellis, N. Fornengo, G. Mignola and S. Scopel, 
"Searching for relic neutralinos using neutrino telescopes," Astropart. Phys. 5, 
333 (1996) [hep-ph/9603342|. 

[89] W. H. Press and D. N. Spergel, "Capture by the sun of a galactic population of 
weakly interacting, massive particles," Astrophys. J. 296, 679 (1985). 

[90] J. Silk, K. Olive and M. Srednicki, "The photino, the sun, and high-energy neu- 
trinos," Phys. Rev. Lett. 55, 257 (1985). 

[91] M. Srednicki, K. A. Olive and J. Silk, "High-Energy Neutrinos From The Sun 
And Cold Dark Matter," Nucl. Phys. B279, 804 (1987). 

[92] J. S. Hagelin, K. W. Ng and K. A. Olive, "A High-Energy Neutrino Signature 
From Supersymmetric Relics," Phys. Lett. B180, 375 (1986). 

[93] K. Ng, K. A. Olive and M. Srednicki, "Dark Matter Induced Neutrinos From The 
Sun: Theory Versus Experiment," Phys. Lett. B188, 138 (1987). 

[94] J. Ellis and R. A. Flores, "Realistic Predictions For The Detection Of Supersym- 
metric Dark Matter," Nucl. Phys. B307, 883 (1988). 

[95] J. L. Feng, K. T. Matchev and F. Wilczek, "Prospects for indirect detection of 
neutralino dark matter," Phys. Rev. D 63, 045024 (2001) [astro-ph/00081 15 1. 

[96] J. D. Lykken and K. T. Matchev, "Supersymmetry signatures with tau jets at the 
Tevatron," Phys. Rev. D61, 015001 (2000) [hep-ph/9903238|; "Tau jet signals 
for supersymmetry at the Tevatron," hep-ex/99 10033, 



62 



[97] K. T. Matchev and D. M. Pierce, "Supersymmetry reach of the Tevatron via 
trilepton, like-sign dilepton and dilepton plus tau jet signatures," Phys. Rev. D60, 
075004(1999) |hep-ph/9 9T3Z52| . 

[98] K. T. Matchev and D. M. Pierce, "New backgrounds in trilepton, dilepton and 
dilepton plus tau jet SUSY signals at the Tevatron," Phys. Lett. B467, 225 (1999) 
[hep-ph/99073U5). 

[99] G. F. Giudice and R. Rattazzi, "Theories with gauge-mediated supersymmetry 
breaking," Phys. Rept. 322, 419 (1999) [hep-ph/980f27T|. 

[100] H. Pagels and J. R. Primack, "Supersymmetry, Cosmology And New Tev 
Physics," Phys. Rev. Lett. 48, 223 (1982). 

[101] S. Weinberg, "Cosmological Constraints On The Scale Of Supersymmetry 
Breaking," Phys. Rev. Lett. 48, 1303 (1982). 

[102] L. M. Krauss, "New Constraints On Tno' Masses From Cosmology. 1. Super- 
symmetric Tnos'," Nucl. Phys. B 227, 556 (1983). 

[103] D. V. Nanopoulos, K. A. Olive and M. Srednicki, "After Primordial Inflation," 
Phys. Lett. B 127, 30 (1983). 

[104] M. Y. Khlopov and A. D. Linde, "Is It Easy To Save The Gravitino?," Phys. Lett. 
B 138 (1984) 265. 

[105] J. R. Ellis, J. E. Kim and D. V. Nanopoulos, "Cosmological Gravitino Regener- 
ation And Decay," Phys. Lett. B 145, 181 (1984). 

[106] R. Juszkiewicz, J. Silk and A. Stebbins, "Constraints On Cosmologically Regen- 
erated Gravitinos," Phys. Lett. B 158, 463 (1985). 

[107] M. Bolz, A. Brandenburg and W. Buchmuller, "Thermal production of graviti- 
nos," Nucl. Phys. B 606, 518 (2001) [hep-ph/0012052|. 

[108] J. L. Feng, A. Rajaraman and F. Takayama, "Superweakly-interacting massive 
particles," Phys. Rev. Lett. 91, 011302 (2003) [hep-ph/0302215|. 

[109] J. L. Feng, A. Rajaraman and F. Takayama, "SuperWIMP dark matter signals 
from the early Universe," Phys. Rev. D 68, 063504 (2003) [hep-ph/0306024|. 

[110] L. Covi, J. E. Kim and L. Roszkowski, "Axinos as cold dark matter," Phys. Rev. 
Lett. 82, 4180 (1999) [hep-ph/99052L2|. 

[Ill] L. Covi, H. B. Kim, J. E. Kim and L. Roszkowski, "Axinos as dark matter," JHEP 
0105, 033 (2001) [hep-ph/0101009]. 



63 



[112] L. Covi, L. Roszkowski, R. Ruiz de Austri and M. Small, "Axino dark matter 
and the CMSSM," hep-ph/ 0402240| 

[113] D. Hooper and L. T. Wang, "Evidence for axino dark matter in the galactic 
bulge," |hep-ph/0402227Jl 

[114] J. L. Feng, A. Rajaraman and F. Takayama, "Graviton cosmology in universal 
extra dimensions," Phys. Rev. D 68, 085018 (2003) |hep-ph/0307375|. 

[115] J. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, "Gravitino dark matter in the 
CMSSM," |hep-ph/03T2262| 

[116] J. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, "Very constrained minimal 
supersymmetric standard models," hep-ph/04051 10 

[117] J. L. Feng and T. Moroi, "Tevatron signatures of long-lived charged sleptons 
in gauge-mediated supersymmetry breaking models," Phys. Rev. D 58, 035001 
(1998) |hep-ph/9712499| |. 

[118] D. Acosta, talk given at the 14th Topical Conference on Hadron Collider Physics, 
September 29 - October 4, 2002, Germany. 

[119] X. Chen and M. Kamionkowski, "Particle decays during the cosmic dark ages," 
|astro-ph/03 10473| 

[120] K. Sigurdson and M. Kamionkowski, "Charged-particle decay and suppression 
of small-scale power," astro-ph/031 1486 

[121] J. R. Ellis, D. V. Nanopoulos and S. Sarkar, "The Cosmology Of Decaying Grav- 
itinos," Nucl. Phys. B 259, 175 (1985). 

[122] J. R. Ellis, G. B. Gelmini, J. L. Lopez, D. V. Nanopoulos and S. Sarkar, "Astro- 
physical Constraints On Massive Unstable Neutral Relic Particles," Nucl. Phys. 
B 373, 399 (1992). 

[123] M. Kawasaki and T. Moroi, "Electromagnetic cascade in the early universe and 
its application to the big bang nucleosynthesis," Astrophys. J. 452, 506 (1995) 
(astro^ph/9412055| |. 

[124] E. Holtmann, M. Kawasaki, K. Kohri and T. Moroi, "Radiative decay of a long- 
lived particle and big -bang nucleosynthesis," Phys. Rev. D 60, 023506 (1999) 
[hep-ph/98055P5|. 

[125] M. Kawasaki, K. Kohri and T. Moroi, "Radiative decay of a massive particle 
and the non-thermal process in primordial nucleosynthesis," Phys. Rev. D 63, 
103502 (2001) [hep-ph/0012279|. 



64 



[126] T. Asaka, J. Hashiba, M. Kawasaki and T. Yanagida, "Spectrum of back- 
ground X-rays from moduli dark matter," Phys. Rev. D 58, 023507 (1998) 
|hep-ph/9802271|. 

[127] R. H. Cyburt, J. Ellis, B. D. Fields and K. A. Olive, "Updated nucleosynthe- 
sis constraints on unstable relic particles," Phys. Rev. D 67, 103521 (2003) 
|astrcTph/0211258| . 

[128] M. H. Reno and D. Seckel, "Primordial Nucleosynthesis: The Effects Of Inject- 
ing Hadrons," Phys. Rev. D 37, 3441 (1988). 

[129] S. Dimopoulos, R. Esmailzadeh, L. J. Hall and G. D. Starkman, "Limits On Late 
Decaying Particles From Nucleosynthesis," Nucl. Phys. B 311, 699 (1989). 

[130] M. Y. Khlopov, Cosmoparticle Physics, Singapore: World Scientific, 1999. 

[131] K. Kohri, "Primordial nucleosynthesis and hadronic decay of a massive par- 
ticle with a relatively short lifetime," Phys. Rev. D 64, 043515 (2001) 
|astro-ph/0103411|. 

[132] K. Hagiwara et al. [Particle Data Group Collaboration], "Review Of Particle 
Physics," Phys. Rev. D 66, 010001 (2002). 

[133] R. H. Cyburt, B. D. Fields and K. A. Olive, "Primordial Nucleosynthesis in Light 
of WMAP," Phys. Lett. B 567, 227 (2003) [astro-ph/0302431 1. 

[134] D. Kirkman, D. Tytler, N. Suzuki, J. M. O'Meara and D. Lubin, "The cosmolog- 
ical baryon density from the deuterium to hydrogen ratio towards QSO absorp- 
tion systems: D/H towards Q1243+3047," Astrophys. J. Suppl. 149, 1 (2003) 
[astro-ph/0302006]. 

[135] S. Buries, K. M. Nollett and M. S. Turner, "Big-Bang Nucleosynthe- 
sis Predictions for Precision Cosmology," Astrophys. J. 552, LI (2001) 
JastrcTph/0010171| . 

[136] K. Jedamzik, "Lithium-6: A Probe of the Early Universe," Phys. Rev. Lett. 84, 
3248 (2000) [astro-ph/9909445]. 

[137] K. Jedamzik, "Did something decay, evaporate, or annihilate during big bang 
nucleosynthesis?," astro-ph/0402344 



[138] M. Kawasaki, K. Kohri and T. Moroi, "Hadronic decay of late-decaying particles 



and big-bang nucleosynthesis," astro-ph/0402490, 



65 



[139] J. L. Feng, S. f. Su and F. Takayama, "SuperWIMP gravitino dark matter from 
slepton and sneutrino decays," hep-ph/0404198, 

[140] J. L. Feng, S. Su and F Takayama, "Supergravity with a gravitino LSP," 
|hep-ph/0404231| 

[141] W. Buchmuller, K. Hamaguchi, M. Ratz and T. Yanagida, "Supergravity at col- 
liders," Phys. Lett. B 588, 90 (2004) [hep-ph/D3U2T73). 

[142] W. Buchmuller, K. Hamaguchi, M. Ratz and T. Yanagida, "Gravitino and gold- 
stino at colliders," hep-ph/0403203 . 

[143] F. Wang and J. M. Yang, "SuperWIMP dark matter scenario in light of WMAP," 
hep-ph/0405186 

[144] W. Hu and J. Silk, "Thermalization Constraints And Spectral Distortions For 
Massive Unstable Relic Particles," Phys. Rev. Lett. 70, 2661 (1993). 

[145] D. J. Fixsen et ah, "The Cosmic Microwave Background Spectrum from the Full 
COBE/FIRAS Data Set," Astrophys. J. 473, 576 (1996) [astro-ph/9605054|. 

[146] http://map.gsfc.nasa.gov/DIMES 

[147] J. L. Feng and M. M. Nojiri, "Supersymmetry and the linear collider," 
|hep-ph/0210390| 

[148] http://www.rssd.esa.int/index.php ?project=PLANCK 

[149] M. Drees, Y. G. Kim, M. M. Nojiri, D. Toya, K. Hasuko and T. Kobayashi, 
"Scrutinizing LSP dark matter at the LHC," Phys. Rev. D 63, 035008 (2001) 
[hep-ph/000 7202| . 

[150] J. L. Feng, M. E. Peskin, H. Murayama and X. Tata, "Testing supersymmetry at 
the next linear collider," Phys. Rev. D 52, 1418 (1995) [hep-ph/9502260]. 

[151] R. D. Peccei and H. R. Quinn, "Constraints Imposed By CP Conservation In The 
Presence Of Instantons," Phys. Rev. D 16, 1791 (1977). 

[152] S. Weinberg, "A New Light Boson?," Phys. Rev. Lett. 40, 223 (1978). 

[153] F. Wilczek, "Problem Of Strong P And T Invariance In The Presence Of Instan- 
tons," Phys. Rev. Lett. 40, 279 (1978). 

[154] A. Kusenko and M. E. Shaposhnikov, "Supersymmetric Q-balls as dark matter," 
Phys. Lett. B 418, 46 (1998) [hep-ph/9709492|. 



66 



[155] A. Kusenko, V. Kuzmin, M. E. Shaposhnikov and P. G. Tinyakov, "Experimental 
signatures of supersymmetric dark-matter Q-balls," Phys. Rev. Lett. 80, 3185 
(1998) |hep-pfa/97122Hl |. 

[156] K. Enqvist and A. Mazumdar, "Cosmological consequences of MSSM flat di- 
rections," Phys. Rept. 380, 99 (2003) [hep-ph/0209244|. 



67 



